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NUMERICAL SIMULATION OF THE WORLD OCEAN CIRCULATION 
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by 

Kenzo Takano 7 Yale Mintz and Young June Han 
University of California/ Los Angeles 

Abstract 

A multi-level model, based on the primitive equations, is developed for 
simulating the temperature and velocity fields produced in the world ocean by 
differential heating and surface wind stress. The model ocean has constant depth, 
free slip at the lower boundary, and neglects momentum advection; so that there 
is no energy exchange between the barotropic and baroclinic components of the 
motion, although the former influences the latter through temperature advection. 

The ocean model was designed to be coupled to the UCLA atmospheric 
general circulation model, for the study of the dynamics of climate and climate 
changes. But here, the model is tested by prescribing the observed seasonally 
varying surface wind stress and the incident solar radiation, the surface air 
temperature and humidity, cloudiness and the surface wind speed, which, together 
with the predicted ocean surface temperature, determine the surface flux of radiant 
energy, sensible heat and latent heat. From an initial state, in which the temper- 
ature as a function of depth and latitude is given, the integration is carried forward 
over several decades of simulated time, which is longer than ihe adjustment iime 
scale for the upper levels of the model. In this test, the model successfully simulates 
most of the large scale features of the world ocean circulation, including a good 
simulation of the ocean surface temperature. 


Permanent affiliation: Institute of Physical and Chemical Research, Wako-shi, 

Saitama-ken, Japan. 
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Introduction 


Governing equations 

With momentum advection neglected, and with the Boussinesq approxi- 
mation for the density and incompressibility, the governing equations are 


lr = - — vp -2flxw+r , 
3t Pq 

{1) 

dp _ 

£-~ 9p ’ 

(2) 

p = P o 0-a|T-T o |), 

(3) 

|i + v.aw) + |j(Tw) = | , 

(4) 


V* W 


+ ?w 

dz 


55 0 , 


( 5 ) 
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where ty is horizontal velocity, w is vertical velocity, p is pressure, p is 
density and p is a constant, T is temperature and is a constant, /? is the 
earth's rotation vector, g is the acceleration of gravity, a is the coefficient 
of thermal expansion, and c is the specific heat of sea water, z is the positive 
upward vertical coordinate and t is time. IF is the horizontal frictional force 
per unit mass and Q is the rate of heating per unit mass, which are taken as 

' ■ (6) 

and 


Q 


= c A|_|V® T + 


cK 9 S T 
6 9z a ' 


(7) 


where A^ is the coefficient of horizontal eddy viscosity, Is the coefficient 
of horizontal eddy heat diffusion, and K is the coefficient of both the vertical 
eddy viscosity and the vertical eddy heat diffusion. The coefficient 6 (in 7) is 
defined as 



when 


st r^o 

9zl<0 , 


as a parameterization of strong vertical mixing that restores a neutral 
stratification whenever a vertically unstable stratification developes. 

The boundary conditions at the ocean bottom (z = H) are 

w = 0 , 



/ 



( 8 ) 

(9) 


and 


( 10 ) 
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so that there is no vertical motion, no momentum transfer and no heat tiansfer 
across the lower boundary. 

The boundary conditions at the ocean surface (z = 0) are 


w = 0 , 


(10 


and 


3 w _ r s 

3z p K ' 
r o 


§1 = 
3z 


Q 

s 

cKp 


/ 


( 12 ) 


(13) 


where is the vector wind stress and Q is the downward heat flux at the 
$ s 

ocean surface. Setting w = 0 at the ocean surface filters out external gravity 
waves and allows the use of a long computational time step. 


The boundary conditions at the lateral walls are zero normal velocity 
and zero tangential velocity, and zero heat transfer through the walls. 


Calculational procedure 

The assumption of no momentum advection, constant depth, and free slip 
at the lower boundary, decouples the barotropic component of the velocity (the 
vertical mean of the horizontal velocity) from the baroclinic component of the 
velocity (the deviation of the horizontal velocity from its vertical mean). The 
former Is governed entirely by the wind stress curl and is independent of the 
heating, while the latter is governed by both the stress and the heating. There 
is no energy exchange between the two components, although the barotropic 
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component affects the baroclinic component through the horizontal temperature 
advection (Eq. 4)* 


The barotropic component of the velocity is obtained by taking the 
curl of the vertical integral of Eq. 1 , or 


where 


and 


3/M 

3t 


l 

- ~ vP -2ftx/M + -i- + A x ,v*/M , 

o p M 

*o 



(14) 


(15) 


where Ikx v<f> = /M defines the stream function <£> of the vertically integrated 
velocity IM, f = 2 n sirup is the coriolis parameter, X is longitude and <p is 
latitude, Ik is unit vertical vector, and a is the Earth's radius. This vorticity 
equation, for the barotropic component of the velocity. Is solved as an initial- 
value boundary -value problem in the way which is described below. 


The baroclinic component of the velocity is obtained by substituting for 
W='W' + W-’ and vp = ’vp 1+ (vp) 1 in Eq. 1, where the bar denotes the 

vertical mean from z = -H to z “ 0, and the prime denotes the deviation from the 
vertical mean. Then, subtracting the vertical mean of Eq. 1 , we obtain 

If = - ' 2flx +a m^ w ‘ + ( k sr- ■ h ) , < 16) 

o o 

which is the prediction equation for the baroclinic component of the velocity. 
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The temperature is obtained from Eq. 4. 

(v p) 1 is found by taking the horizontal gradient of Eq. 2, 3(vp)/3z = .-gVp, 
integrating over the entire depth of the ocean from an arbitrary value of vp 
at z — -H , and subtracting the vertical mean, / vp\ This insures that the vertical 
mean of (vp) 1 is identically zero. 

The vertical velocity w is calculated by Eq. 5. 

Grid structure, domain, and finite difference scheme 

The finite difference calculations are made using spherical coordinates 
and the same grid structure that was used by Haney (1971a). As shown in Figure 1 , 
the T, p, <$> , and w points are located at the coast lines. The W and W* 
points match the distribution of the horizontal velocity points 
of the UCLA atmospheric general circulation model, and therefore include points 
on the equator. In the atmospheric model the horizontal grid size is 4 in latitude 
and 5° in longitude; but in the ocean model the zonal grid size is reduced to 2.5° 
in longitude to obtain a somewhat better resolution of the western boundary 
' currents. 


The model ocean is divided into five layers, which lie between 0, 70, 

380 , 960 , 2020 and 4000 m in depth, and the horizontal velocities and temper- 
atures are predicted for each layer. To calculate the vertical eddy diffusions of 
momentum and heat, these velocities and temperatures are attributed to the 20, 120, 
640, 1280 and 2760 m depths. The vertical velocities are calculated at the interfaces 
of the five layers. 

The horizontal domain extends from 66°S to 74°N, which excludes the 
Arctic Ocean. The northern boundary follows the edge of the continental shelf of 
the composite landmass of North and South America, Europe, Africa and Asia. 
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The southern boundary approximately follows the edge of the continental shelf 
of Antarctica . New Zealand, and Australia combined with New Guinea, are 
retained as islands. 


The finite differencing scheme is the same as the one used by Haney 
(1971a), except that his nine point specification of the vorticity of the barotropic 
velocity component, which was almost the same as a 45° rotated five point scheme, 
is here replaced by the usual nearest five point scheme. This prevents a com- 
putationally false two grid point noise from being produced by the irregular lateral 
boundaries and irregular wind stress. 

The time integration is carried forward by the leap-frog scheme, with 
the Euler-backward scheme applied every five time steps* 

Prescribed parameters 

The prescribed constant parameters are given in Table 1. 


Table 1 . Prescribed Constants 


2 

g = 980 cm/sec 
0 ” 7.29 x 10 sec 

o 

R = 6.37 x 10 cm 
c = 0.93 cal/cm/degr 
a = 2.5 x 10 degr 

— lgm/cm^ 


T = 
o 

H - 



K = 
At = 


0°C 

4 x 10 5 cm 

in 9 2 / 

ID cm /sec 

2.5 x 10^cm^/sec 

1 cm^/sec 

40 min, 8 hours 


40 minutes is the time step for the integration of the vorticity equation for the 
barotropic velocity component, and 8 hours is the time step for the integration 
of the primitive equations for the baroclinic velocity component. 
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When the oceanic model will be coupled to the atmospheric general 
circulation model/ the surface wind stress will be determined by the atmospheric 
model. But/ here, for the purpose of testing the oceanic model/ the surface 
wind stress is taken as the observed seasonally varying stress given by Helierman 
(1967/ 1968). Hellerman's values of the mean wind stress for the four seasons of 
the year (Dec ~Jan -Feb; Mar-Apr-May; June-July-Aug; Sept "Oct- Nov) are 
subjectively interpolated and extrapolated to cover the domain of the oceanic 
model. A 1-2-1 zonal filter is used to remove some of the two grid interval 
noise in Hellerman's tabulations. Then, at each grid point, the four seasonal 
values of ^ are used to calculate the annual mean and the amplitude and phase 
of the first annual harmonic of the stress; and from these parameters the surface 
wind stress is prescribed as a time varying function of the model year. 

Q , the net downward heat flux at the oc ean surface, has four components, 
s 

Q ~ Q - (Q+G +Q.) , (17) 

s I B F L 

where Q is the net downward flux of solar radiation, Q D is the net upward flux 
I o 

of infrared radiation, Qp is the net upward flux of sensible heat, and is the 
net upward flux of latent heat. 

When the oceanic model will be coupled to the atmospheric general 
circulation model, the net fluxes of solar and infrared radiation at the ocean 
surface will depend on the atmospheric state parameters in the fairly complex 
way given by Katayama (1973). In the present test of the ocean model, 
however, we use the relatively simple formulations, given by Angstrom (1922) 
for the solar radiation flux. 


Q I =Q I (1 -0.7|)(1 -A s ) , 
o 


( 


( 18 ) 
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and by Brunt (1952) for the infrared radiation flux, 

Q =oT 4 x 0.985 [0.39 - 0.05 e^]x [1 -0.6 (C/8) 9 ] 

D S 

Where Q is the direct plus indirect solar radiation reaching the ground through 
■*o 

a cloudless atmosphere, C is the number of octants of the sky covered by clouds, 
Ag is the albedo of the ocean surface, cr is the Stefan-Boltzman constant; and 
e is the vapor pressure of the air (in mb) at 10 m above the ocean surface. 


09 ) 


The annual mean and the amplitude and phase of the first annual harmonic 

of Q are obtained from the Smithsonian Meteorological Tables (1966, fables 133 
lo 

and 136), with an atmospheric transmission coefficient of 0.7. The albedo of the 
ocean surface is taken as 0.07. 

The annual mean and the amplitude and phase of the first annua! harmonic 
of the number of octants of the sky covered by clouds, C, are calculated from the 
charts in the cloud atlas of Miller (1971). 

The annual mean and the amplitude and phase of the first annual harmonic 
of the atmospheric vapor pressure, e, are calculated from tabulations of the surface 
air dewpoint temperature in the atlases of Taljaard et. al. (1969) and Crutcher and 
Meserve (1970). 

These time varying parameters are used, in Eqs. 18 and 19, to prescribe 
the net flux of solar radiation and to calculate the net flux of infrared radiation 
at the ocean surface . 

For the sensible and latent heat fluxes, the formulations for the test 
calculations are nearly the same as those that will be used with the combined 
atmospheric and oceanic models, 

Q F = p A C D C pl \l ( V V ' ( 2 °) 


and 


Q L = PA C D L I^Ai (q *' q A ) ' 


(21) 
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where is the surface air density/ C^(= 1 *23 x 10 ) is the surface drag 

coefficient, C is the specific heat of air, | W A | is the surface air speed, L 
P * . ^ . . 

is the latent heat of e vaporization, q^ is the saturation water vapor mixing 

ratio at the ocean surface temperature, and q^ is the water vapor mixing ratio 

of the surface air. 

The annual means and the phases and amplitudes of the first annual 
harmonics of and q^ are calculated from tabulations of the surface air 
temperature and the surface air dewpoint temperature in the atlases of Taliaard 
et. al. (1969) and Crutcher and Meserve (1970). The annual mean and the 
phase and amplitude of the first annual harmonic of the surface (scalar) wind 
speed are calculated from data in the atlas of MacDonald (1934, charts 27 
through 30). ^ 


The simulated barotropic circulation 


The simulated baroclinic circulation 



The simulated surface heat flux 


Summary and conclusions 
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BAROTROPIC CIRCULATION IN A MULT I- CONNECT ED GLOBAL OCEAN 

fipp=0 IF THE FOUR NE ARES R SURROUNDING POINTS A R.E WATER POINTS, 

= 1 if ONE OF THE FOUR NEAREST SURROUNDING POINTS IS ON THE 
‘ COAST OF T HE WAIN LANDMSS (AFRICA, EUROPE, ASIA), 

= 3 IN SIDE THE ANTARCTICA, 

=30 ON THE COAST OF THE ANTARCTICA, 

=4 INSIDE AUSTRALIA, 

=40 ON THE COAST OF AUSTRALIA, 

=5 INSIDE N EH ZEALAND, 

=50 ON THE COAST OF N5WZEAL AND. 

MULTIPLY THE STREAM FUNCTION VALUES BY R*H ( B= EARTH 1 S 
RADIUS, . H= OCEAN DEPTH) TO GET THE STREAM FUNCTION VALUES 
FOR THE VOLUME TRANSPORT (INTEGRAL OF THE HORIZONTAL VELOCITY 
OVER THE W HOLE DEPTH. ; 

HQQ ~ 1 FOR THE EATER POINTS, 

=2 FOR THE MAIN LAND MASS, 

=3 FOR THE ANTARCTICA, 

=4 FOR AUSTRALIA, 

= 5 FOR NENZB ALA ND . 

THE SUBROUTINE HAJII4E SPECIFY THE SHAPE OF THE OCEAN (MPP AND 
MQQ) , THE KIND STRESS TO COMPUTE ITS CURL AND OTHER CONSTANTS 
AS NELL AS THE INITIAL FIELD OF THE STREAM FUNCTION. 

THE SUBROUTINE RESID GIVES TEE CHANGE OF TEE SUM OF THE 
RESIDUALS (RES) AT EACH SEPARATE LAND MASS WHICH TAKES PLACE 
BY THE UNIT INCREASE OF THE STREAM FUNCTION AT EACH SEAPARATE 
LAND MASS POINT, JUST AS IN CASE OF THE BLOCK RELAXATION. 

F AND G DENOTE THE STREAM FUNCTION AT TIME STEP N + 1, N , 
RESPECTIVELY, XU,XV,U, AND V THE BAROTROPIC COMPONENTS AT TIME 
STEPS N + 1 , N, RESPECTIVELY. 

WIND STRESS CD RL=-TW/ ( 2. 0 * Y*X { J) ) / R*H 
X, XM,Y=GRID SIZE 

AH=CO EFFICIENT OF EDDY VISCOSITY 
DT IS THE TIME STEP. 

H IS THE OCEAN DEPTH. CD IS THE TIME (DAY). 

E IS THE WORKING SPACE FOR RELAXATION. 

Ill DENOTES THE EASTERNMOST POINTS TO BE SWEPT. (EX. 146) 

J M DENOTES THE NORTHERNMOST FONTS TO BE SWEPT. 

THE BACKWARD DIFFERENCING SCHEME IN TIME APPLIES EVERY MAT 
TIME STEPS. MATS CONTROLS THE SWITCHING FROM THE LEAP-FROG TO 
BACKWARD OF FROM THE BACKWARD TO LEAP-FROG SCHEMES. 

PRINTING IS DONE EVERY NPRINT TIME STEPS. 

WRINTING IN MAGNETIC TAPE IS DONE EVERY NT APE TIME STEPS. 

WHEN (IN) EQUALS ONE, THE IMAG-TH RECORD IN THE MAGNETIC TAPE 
IS PICKED UP, AND THE INTEGRATION GOES ON UNTIL N N=N END. 

REL= RELAX AT ION ACCURACY 

WIND STRESS CO RL=-TW (I , J) *H (DEPTH ) / (2 . 0*Y*X (J) *RAD) RAD= 
RADIUS OF THE EARTH) 

I N D (1) -30 ANTARCTIC COAST 
IN D ( 2) =4 0 AUSTRALIA COAST 
I ND (3 ) =5 0 NEWZSALAND COAST 
INDIN (1) =3 ANTARCTIC INLAND 
INDIN (2) =4 AUSTRALIA INLAND 
INDIN (3)=5 NEWZEA LA ND INLAND 

FIGURE 1 EXPLAINS THE MEANING OF IX1,LX2, IY 1 AND LY2. 
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FIG. 7 

L,KA,KB SPECIFY THE ST REA 8 FUNCTION POINTS TO BE SWEPT AND 
OTHER LAND POINTS TO BE SKIPPED. 

FIGURE 2 EXPLAINS THEIR MEANINGS » 

* POINT TO BE SWEPT 


J=4 KA(1, J) ***KB(1,J) KA (2, J) *****K3(2,J) K A (3 , J) ****KB ( .1, J) 

L (J) =3 


C FIG. 2 

C 

c 

C THE COMPUTATION STOPS IF TEE NF.RGY EXCEEDS A CERTAIN VALUE, OR 

C IF THE NUMBER OF ITERATION IN RELAXATION EXCEEDS A CERTAIN ? ALU 

C THESE VALUES ARE SPECIFIED IN THE PROGRAM . THEY SHOULD BE CHANG 

C IF NECESSARY. 

C THE COEFFICIENT OF THE OV EE- RE LAX A TION IS GIVEN AT THE BEGIN- 

C NING. 

C 

c 

COMMON AH,DT f Y,H, REL, XM (39) ,FJU(3 9) , X (39) , F ( 1 47 , 39) ,E( 147,3 9) , 

1U {1 47,3 9) ,V(1 47,3 9) ,5 {147, 39) ,XU(147, 39) , XV (147,39) ,TW (147, 39) ,CD, 

2 BET (3 9) , YU ,CC Y, AH , ARY ,RE5 (3) , DDT, 

3 M AT , L (3 9) , LV (39 ) , LA <1 0,37) , LB <1 0, 37) , KA ( 1 0, 37) , KB ( 1 0, 37) , MATS, 

4MPP ( 1 47 , 3 9) , MQQ ( 1 47 ,39) , INDIN (3) ,IND(3) ,NRE5(3) , 

5LX1 (3) , LX 2 (3) ,LY1 (3) ,LY2(3) r If5, 1, JMP1 ,IMM1 

SOB =1.77 5 
SOR=1.B75 
SOR=1.825 
SOR = 1 .85 
SOR 1=SOR- 1.0 
READ (5, 10 00) IN 

1 000 FORMAT (IS) 

IF(IN.EQ. 1) GOTO 1001 
CALL HA J I ME 

CALL B.ESIB 
CD=0. 0 
NN=0 
M AG = 0 

GO TO S55S 

1001 CONTINUE 

READ (5, 1000) I MAG 
IF (IMAG. EQ. 1) GO TO 77 
DO 7 3=2, IMAG 
READ { 10) 

7 CONTINUE 

'7 READ (10) XU,XV,F,G,U,V,TW,HPP,HQQ,L, KA, KB, L V ,LA , LB, AH , DT , H AT, 
1Y,RF:L,XM, FJU, X, BET, I X 1 , LX 2 , LY 1 , LY2 , IN D, INDIN, NN, NEED, RES, 

2NRES, ARY, AR , Y4,E,CCY, MATS , CD, IH , JM, IH El , JMP1 , IMH1, MAG 
5555 CONTINUE 
AH= 1 . OE 9 
AH=3,0K8 

AR=AH/6. 37E8/6.37E8 
ARY= AR/Y 

t n * p ~«m i 


w w 



JiNN = 0„0 

HBAD (5,5) MEND, MAT, NPRINT , NT APE, DT, RED 
^ 5 F0RMAT(4I5,2F,15.7) 

' , WRITE (6,6) I HAG, N END, DT ,R EL, AH, MAT, NT APE, NPRINT 

tf 6 FORMAT (1H ,16HTIME STEP FROM l5 f 2X,3HTO 15, 2 X, 3HDT=E1 5 * 7 , 2 X, 4HREI 
1'EI 5. If 2 X , 3 HA H=E1 5* 7, 2X ,4HMAT=I 3/1 K , 5HTAPE«I4, 5X„ 6HPRI NT=I 4) 



0 DAY=DT/0. 864E5 
*QQ (49,5) - 1 
HQQ(50,5) =1 
MQQ (49,3) =3 
HQQ{50,3)*3 


39 NN=NN+1 


CD=CD+DAY 


NNN = (NN/MAT) * HAT 

I F ( N NN. EC. NN) GO TO 25 

H ATS = 0 

DDT=DT 

DO 1 3=3, OH 

DO 2 1= 1 , IMP1 

R=G (1,3) 

G (I,J)=F(I, J) 

F(I,J)=B 
2 CONTINUE 
1 CONTINUE 
GO TO 26 

25 MATS=1 
PDT=DT*0 - 5 
DO 4 2 J = 3 , 0 M 
DO 43 1=1, I M Pi 
U (I ,J) = XU (I, J) 

V (I, J}= XV (I , J) 

G(I,J)=F (1, J) 

43 CONTINUE 

42 CONTINUE 

26 CONTINUE 

DO 3 0 <7=3, OH 
JP1=J+1 

" JM 1 = J-1 
JM2=J-2 
X«0=XH(0) 

X M 1 = X M ( J M 1 ) 

Cl = Y/XMG 
S 1= Y/XH 1 
CB= X (3P1) /Y 
CS=X (0) /Y 
B = X ( J Hi ) /Y 
BET J = BET (J) *CCY 
LJ=L (J) 

DO 3 1 Ll= 1 ,L J 
11 = KA (LI.,0) 

MM= KB (LL , 0) 

DO 310 I =TI , M M 
IP1=I+1 
IH 1 = 1-1 
I M2=I-2 
UIJ = U (1,0) 

VIJ = V (1,0) 

DPI J-U (IP1, J) 
DH1J-U(IM1,J) 

UIP1 = U(I/0P1) 

UM1 P 1 -U (I Hi , J PI ) 

UIH 1 = U (I , OH 1) 

UH1M1-U (IM1,»7M1) 
UP7H1-U (IP1 ,JH1) 

UIH 2 = U (1 , 0H2) 

V T H2=V(I, JH2) 


33 | 


ft. 



UiiZ J~U J) 

V M2 J = V (IH2, J) 

%. UH1M2=U 
" „ VM1 H 2 = V (T N 1 , J M2) 

* UM2M 1-0 (IM2, JM1) 

VM2M1=V (IM2, JM1) 

VP1 J=V (IP 1 , J) 

Vfll J=V ( I M 1 , J } 

VIP 1 — V (I , JP 1) 

vitn=v (i , jmi) 

VH1P1-V(im # jpi) 

VM1M 1=V ( I Ml, JM1) 

VP1 H 1 =V (IP1 , JM1) 

MP=MPP (I f J) 

IF(MP.EQ. 0) GO TO 304 

IF (MP.EQ. 30) GO TO 304 

IF ( MP *EQ. 50) GO TO 304 

IF { MQQ { IP 1, J) . NE„2) GO TO 27 

UP1 J=-UIJ 

VP1 J=-VI J 

27 IF (MQQ (I r JP1) .NE.2) GO TO 271 
UIP 1=-UIJ 
VIP1=-VI J 

271 IF(MQQ(IMl,JPl) -NE.2) SO TO 272 
UK 1 P 1=-UM 1J 


VM1P1=-VM1J 

272 IF (HQQ(IM2,J) .NS. 2) GO 70 28 
UH2J*-0M 1 J 
V M2 «]--V M1J 

2 8 IF (MQQ (I M2, JM 1) .NE. 2) SO TO 2 9 
UM2M 1=-U B 1M1 
VM2M1-—VM1M1 

29 IF(MQQ(IM1, JM2) -HE, 2) GO TO 291 
UM1W2*-UI1 1H 1 
V Ml H2*-VM1 Ml 

291 IF(KQQ(I, JM2) • NE. 2) GO TO 292 
UIB2=-UIM1 

VIM2=-VIM1 

292 IF (MQQ (IP 1, JM1) .NE. 2) GO TO 303 
U P 1 M 1=-DIM1 

VP1 M1~- VI Ml 


303 CONTINUE 

3 0 4 DF=— BET J * (XV (I, J) +XV (IM1 ,J) +XV +XV (I, JM1) } *0.25 

5S= (Cl* (0 31 J+OIJ-UM2J-U Pi J) - (UI Pi — U IJ +0 M IP 1— U M 1 J) *CB+2. 0*CS* 

1 (UIJ-UIM 1 + UM1 J-IJM1M 1) *S1* (UP1M1 +UM2M1 - UIM1-OM 1M1) - (UIfU-UIM2+TJf! 1M1 
2-UMlM2)*B ) *ARY4TW (I, J) +DF 
DF=CS*( VI J-VI M1-V M1J+VM1M 1) 

SS= AR/XMJ* {Cl* (-3. 0* (VI J-VM1 J) + V P 1 J-VM2 J ) ~DF 
1 + (V IP 1-VIJ- VM 1P1 + VM1J) *CB) +AR/XM1 * (Si * (3. 0* ( VMM 1- VIM 1) +VP1 Ml 
2-VM2M1) + (VH1M 1-VH 1H2-VIM1 + VIH2) *B +DF)+SS 
4 E (I , J) =-S S*DDT 
310 CONTINUE 
31 CONTINUE 

E (2, J)=E (IB, J) 

E (IMP1 , J) =E (3 , J) 

30 CONTINUE 


IP=0 



7 IP=IP+1 
**I SL A NDS 
MAP = 0 


C*****ANT ARCTIC A 
C *****MEW2E ALA ND 
C OSTR ALIA 

DO 5021 K-1 , 3 
RES K-RES ( K) 
MP*IH.D (K) 



MpMfc»='iN DIN In) 

II=LX1 (K) 

IG=IX2(K) 

* , IM2=LY1 (K) 

# CH2 = LY2(K) 

R=0 - 0 

DO 5016 J=TM2 f JM2 

DO 5015 1=11 , IG 

IF{HPP (I, J) - NE. MP) GO TO 5015 

E=H + E (I, J) 

5015 CONTINUE 

5016 CONTINUE 

IF(ABS( H ) . L E. RE L) GO TO 5021 
BF=-IVRESK 
HAR=«AR+ 1 
IF(R.EQ.I) IM2=2 
DO 5022 J=IH2 , JM2 
JP1=0+1 

CS= Y/X ( J) 

C B= 2.0* (CS + 1 . 0/CS) 

S 5= X (OP 1 ) /Y 
C1 = X (JM 1) /Y 
DO 5023 I =11, IG 
MPPIJ=MPP (1,0) 

HMMM= (K PPIO-MP) * (8PPIJ-HPMP) 

XF ( MKMK. NE. 0) GO TO 502 3 
F (I, J)=F (I, J) + DE 
I ill = 1-1 
IP1 =1 + 1 
CCC=DF*CS 

E (I , JP1) =E (I, JP1) +SS*DF 
E {I, JM1 ) =E (I, JM1) +C 1 * DF 
£ (IP 1 , J ) = E (TP1,J) +CCC 
E (IK 1 r J) =E (IK 1, J) +CCC 
E (I, 0)=E (1,0) -CB*DF 

5023 CONTINUE 
5022 CONTINUE 

IF(K.EQ.2) GO TO 5021 
IF(K.EQ.I) GO TO 5024 
DO 5026 J=8, 1 2 
E (2 , J) = E(IM,J) 

5026 CONTINUE 

GO TO 5021 

5024 CONTINUE 
E (IMP1, 5) =E (3,5) 

E (2, 5)=E(TH,5) 

CCC=DF*Y/X (4) 

E (3 , 4 ) =E (3,4) +CCC 
E (IN, 4) =E (IM, 4) +CCC 
E (I M P 1 , 4) =E (3,4) 

E (2 , 4) -E ( IM , 4 ) 

5021 CONTINUE 
C*****NEWZE ALAND 
C*****ANTARCTICA 

DO 32 J=4,OM 
OM1 =0-1 
JPl=J+1 
X0=X{0) 

CS= Y/XJ 
SS= X (OP1) /Y 

c 1 = x ( JM 1 ) / Y 

CB=2. 0* (CS + XJ/Y) 

LO=L (J) 

DO 33 LL= 1 , LO 
T T= K A (LL , J) 



uu-cvJJ ll»L, J) 

DO 330 I *11, KM 
IF(MPP(I r J) «GE.3) GO TO 330 
r »=E(I,J) 

IF (ABS(B) . LE. .PEI) GO TO 330 

KAR=MAR+1 

D?= F/CB*SOR 

E (T ,0) =-R*SOS 1 

F (I, J)=F (I, J) +DF 

IPl = I + 1 

CC C=DF*CS 

E (IP 1 ,3) = E (IP 1 , J) +CCC 
I Ml =1-1 

E (I HI ,J) =E(IK1, J) +CCC 
E (I, JP1 ) =E (I, JP1) +SS*DF 
E (I, JM1) =E (I, JK1) + C1*DF 
IF (I.NE.3) GO TO 331 
E {I HP1, JP1) =E (3, JP1) 

E (IMP1, JM1) = E (3, JM1) 

E{IMP1, J) =E (3,J) 

E(IM,J) =2 (2,J) 

GO TO 330 

331 IF { I . N E„ IM) GO TO 332 
E<2, J)*E(IM,J) 

E (2, JP1) =E (IM ,JP1) 

E {2, JM1 )*E (IM,JM1) 

E {3, J)*35 (IMP1 , J) 

GO TO 330 

332 IF ( I - NE . 4 ) GO TO 333 
E (IH P 1 , J ) =E (3, J) 

GO TO 330 

333 IF (I.NE.IMM1) GO TO 330 
E (2, J)=E (IM,J) 

330 CONTINUE 
33 CONTINUE 

32 CONTINUE 

IF (MAE. EQ. 0} GO TO 36 

IF (IP. GT. 1000) GO TO 1112 

IF { IP. RE. 1) GO TO 37 
WRITE (6, 38) IP, MAS 

33 F OE M A T ( 1 H , 7 X ,6 HKAISU=I 5, 5 X, 1 2HTEN NO KAZU=X5) 

GO TO 37 

36 EN = 0 . C 

WBITB(6,3fl) IP, MAR. 

DO 70 J*3,JK 
F (1, J) =F (III L1 1 ,J) 

F(2,J)=F(IM,J) 

F(IHP1, J) =F (3,J) 

70 CONTINUE 

DO 23 J*3,JM 
JM1=J-1 
Jp1=J+1 
XMJ=XM (J) * 2. 0 
DF=XMJ*Y4 
I,VJ = LV (J) 

DO 24 L J = 1 , LVJ 
II - LA (LJ,J) 

IG=LB (LJ, J) 

DO 244 I=II, I G 
IP1=I + 1 

0 (I, J)=XU (I, J) 

V (I,J)=XV (I,J) 

XU (I ,0) = 

1 <F (I, J)+F (IP1, J) -F (I, JP1 ) - F(IP 1, JP1) ) /CC Y 

XV(I,J) = 

1 (F (IP1 , J) +F (IP1, JP 1) -F (I ,J)-F (I, JP1) )/XMJ 



AiW* tSN+JJr^ ( XU (1, d) --^2} 

244 CONTINUE 
, 24 CONTINUE 
♦ XU ( I MP1 , 3 ) -XU ( 3 , 3 ) 

XV(IMP1, J)=XV (3,J) 
l) (IWP1, J) =U{3,J) 

V(IMP1, J) = V ( 3 , 3 ) 

XU (2 ,3) =XU (114,1) 

xv( 2 , 0 )=xv (in r J) 

XU( 1 , 3 ) = XU (IMM1, J) 

XV (1 , 3 ) = XV (IMM1 , «J) 

U (2, J)=U (IM,0) 

V (2 , J) =V ( IH , J ) 

IT (1 r J) = U ( I K M 1 , 3 ) 

V (1 , J)=V (IMK1 , 3 ) 

23 CONTINUE 

WRITE (6,34) NN,CD,EN,F (1,3) 

34 FOR MAT ( 1 H , 5HSTEP-I5, 5X , 5 HTI ME=E1 5. 7, 5X, 7HENESGY=E 1 5. 7, 5X,4HAC 
1.4) 

IF (EN.GT-5.OE1) GO TO 1111 

IF (HU.GE. NEND) GO TO 88 

IF (MATS. I? E. 0) GO TO 10 

IF (NN.NE. ( (NN/NTAPE) * NT. APE) ) GO TO 9 

MAG=HAG+1 

WRITE (10) XU,XV,F ,G,U, V,TW, MPP,MQQ,L, KA,KB,LV,lA,I,B r AH,DT,MAT, 
1 Y ,REI ,XM,FJU, X,BET,LX1 , LX 2 , L Yl, 1Y 2, IN D, IN DIN, NN, N END, RES, 

2NBES , AR Y , AS , Y4 , E , CC Y, H ATS , CD,IM, J M, IM Pi # JMP 1 # IHM1 # M AG 
WRITE (6, 8) HAG 

9 IF (NN .EQ. ( (NN/HPHINT) *NPRINT) ) GO TO 88 
10 ENN=EN 

IF (H ATS. EQ. 0) GO TO 3 9 
M AT S— HATS+ 1 

IF (MATS. EQ. 3) GO TO 4 6 
DO 44 J=3,JM 
DO 4 5 3 — 1 , IMP 1 
F(I,J)«G (1,0) 

U (I , J} = XU (1,1) 

" V (I, J)=xv (I, J) 

45 CONTINUE 
44 CONTINUE 

GO TO 26 

46 DO 47 J=3,JM 
DO 48 1=1, IMP 1 
R=F (I,J) 

F(I,J) = G (I,J) 

G (I , J) -R 
48 CONTINUE 

47 CONTINUE 

DO 49 1=3,0 M 
JP1=J*1 
XMJ-XK(J) *2.0 
L VJ-IV (0) 

DO 55 L J = 1, L V J 
I I = L A (L J , 0) 

IG=LB <IJ, 3 ) 

DO 550 I -II , I G 
IP1=I + 1 
I U(I,J) = 

1 (F (I, J) +F (IP1, J) -F (I, JP1 )-F (IP 1, JP 1) ) /CCY 

V (I , 3 ) = 

1 (F (IP1 , J) + F{IP1, JP1 ) -F (I,J) “F (I, JP1) ) /XMJ 

550 CONTINUE 
55 CONTINUE 
U (1,J)=U 

V {1, J)=V <IMH1,J) 
n (?, J) = U (TM,J) 


3 if 


V(2 f J)“VU tt,J) 

0 { I H P 1 , J ) = U (3,J) 

% V (IBP1, J)=V (3,3) 

v *49 CONTINUE 
‘ f MATS»0 

NN-NN+1 
DDT=DT 
NNT=NNT+ 1 
CD=C D+DA Y 
GO TO 39 
8 8 K=1 
MP=-1 

613 MP=MP+1 
MM-K-MP 
IG-1 8*KP+ 3 

IP(MH.EQ.O) GO TO 614 

II=IG*17 

GO TO 615 

614 1.1=3 8 

615 WRITE (6 ,701) 

701 FORM AT (1 HO, 1 5HSTREAM FUNCTION), 

DO 7 02 J=1,IMP1 

WRIT E (6 , 7 0 3) 0 , (F ( J, I ) , I=IG, II) 

703 FOF.KAT(1I] ,I3,1X, 18F7.4) 

702 CONTINUE 

IF ( MM.NE. 0) GO TO 613 
IF ( NN.LT* NEND) GO TO 10 
M AG=M AG + 1 

WRITE (10) XO,XV,F,G,U,V,TW,MPP,MQQ,L r KA, K 8, L V , L A , LE, AH, DT, MAT, 

1 Y f REL,XM,FJD, X, BET ,LX1,LX2,LY1,IY2, IN D, INDIN, NH,NEND,FES, 

2NRES , ARY, AR ,Y4,£,CCY, BATS ,CD,IM , JM,IMP1 , JHP1 , IHM1 , MAG 
WRITE (6, 3) MAG 

8 FORMAT (1H0,32HN0. OF RECORDS IN MAGNETIC TAPE=I3) 

GO TO 1111 

1112 WRITE (6 , 38) IP, MAE 
1111 STOP 
END 

SUBROUTINE FESID 

RETURN 

END 

SUBROUTINE HAJIME 

RETURN 

END 

SUBROUTINE RE SID 

COMMON AH,DT r Y, H,REL,XM (3 9) ,FJU (3 9) , X (3 9) , F ( 1 47 , 39) , E { 1 47, 3 9) , 

113 (1 47,3 9) , V( 1 47,3 9) ,G ( 1 47,39) r XU ( 14 7, 39) , XV (1 47, 39) ,TN (147,39) , CD, 
2BET (39) ,Y4 ,CC Y, AR , ARY, RES (3) , DDT, 

3 MAT, L (3 9 ) , LV (39) , LA (10,37) , T,B (10,37) , K A ( 1 0, 37 ) , KB ( 1 0 , 37) ,MATS, 

4MPP ( 147, 39) ,KQQ (147,39) , INDIN (3) ,IND< 3) , NRES (3) , 

5LX1 (3), LX 2 (3) , LY 1 ( 3) , LY 2 ( 3) f IM,J8,IMP1,JMPl,IBMl 
DO 1 K= 1 , 3 
DO 8 J~ 1 , JMPl 
DO 9 I=1,IMP1 
E (I , J) = 0* 0 

9 CONTINUE 
° CONTINUE 

s=0.0 
11=1X1 (K) 

1 2=1. X 2 (K ) 

J1=LY1(K) 

J2=IY2 (K) 

MIN=IND(K) 

MINI N=I NDIN (K) 

N U M = 0 

IF(K.EQ.I) J1=2 
DO 3 



TOO 


JP1=J+1 

cs=y/x(j) 
ci=x ( J M 1 5 /Y 
ss = X (JP1)/Y 
CP=-2.0* (CS + 1 .0/CS) 

Do 4 1=11,12 
MP=HPP (I , J) 

MMM= (MP-MIN)* (HP-MININ) 

JF{ MMM.NE, 0) GO TO 4 
IP1 = I + 1 
IM 1 =1-1 
NUM=NUH+1 
E (I, J) = E (I,J) +CB 
E (IP 1 , J) = E (IP 1, J) +CS 
E (I Ml , J) = E{IH1 , J) +CS 
E (I, JP1) =E (I, JP1) +SS 
E(I, JM1) =E {I, Jill) +C1 

4 CONTI NOE 
3 CONTINUE 

IF (K . G£. 2) GO TO 1 0 
J 1=3 

DO 11 J=3 , 3 
CS=Y/X (J) 

E (3, J)=E (3, J) +CS 
E(IU,J)=S (IM, JJ+CS 
11 CONTINUE 
10 CONTINUE 

DO 5 J= J 1 , J2 
DO 6 1=11,12 

IF(MPP(I, J) •EQ.HIN) S=S+E (I , J) 

6 CONTINUE 

5 CONTINUE 
REES (K) =NOM 
RES ( K ) = S 

WRITE (6,7) K,ROM,S 

7 FORE AT ( 1H , 2HK=I2 , 2X, 1 2 HGRID POIH TS=I 3, 2X , 4HRES = E 1 5 . 7) 
"1 CONTINUE 

DO 16 J=1 , 18 
DO 17 1=1, IMP 1 
IF(HQQ(I,J).EQ.1) GO TO 17 
MP'-'MPP (I, J) 

IF (HP.EQ. 2) GO TO 17 
IF (MP.EQ- ( (MP/3) *3) ) GOTO 18 
IF (KP. EQ. ( (MP/4) * 4 ) ) GO TO 19 
HQQ (I,J) =5 
GO TO 17 

18 HQQ (I,J) =3 
GO TO 17 

19 MQQ(I,J) =4 
17 CONTINUE 
16 CONTINUE 

WRITE (6 ,15) (1,1=1,38) 

15 FOR M AT ( 1 H ,4X,38I3) 

DO 12 1=1, IMP 1 

WRITE (6,14) I,(MPF(I,J) ,J=1,38) 

WRIT E (6 , 1 3 ) I, (HQC (I, J) , J=1,38) 

|2 CONTINUE 

*3 FOR M AT ( 1 H ,13, 3X, 3813) 

14 FOR M AT { 1 H ,I3,2X,38I3) 

RETURN 

END 

SUBROUTINE HAJIME 
DATA: 

I H , JM , I M H , JM H 



c\ n o n 



OCfiAN SH Al J ti (SUbKUUIXN r. anAFr, U 

X,DX,D¥ (DX AMD DY ARE THE GRID SIZES IN LONGITUDE AND 
LATITUDE (DEGREES), RESPECTIVELY- Y IS THE LATITUDE 
CORRESPONDING TO 0=0) - 
AH, H , RE L, DT , H AT 

IT* IT- 1 WIND STRESS IS READ TO COMPUTE THE HIND STRESS CURL. 
| * IT=2 HIND STRESS CURL IS GIVEN FROM CARDS, A MAGNETIC 

W TAPE, OR THE SUBROUTINE CURL. 

COMMON AH,DT,Y,H, REL,XM (3 9) , FJU (3 9) ,X (39) ,F(147,39) ,E(147,39) , 

ID (147,39) ,V (1 47,3 9) # 5 (147,39) , XU ( 147, 39) , XV (147,39) ,T.W (147,39) , CD, 

2BET (3 9) , Y4,CCY,AR,ARY,FES (3) , DDT, 

3 MAT,!. (39) ,LV (39) , LA (1 0 , 37) , LB ( 1 0, 37) , KA ( 1 0, 37) , KB ( 1 0,37) ,MATS, 

4MPP( 147, 3 9) ,MQ0<1 47,39) ,INDIN(3) ,IND(3) , USES (3) , 

5 LX 1 ( 3) , LX 2 ( 3) , LY 1 (3) ,LY2(3) ,IH, JH # IKP1, JHP1 r IMHl 
DIMENSION LI (9) , WE (147,3 9) ,WN (147 , 39) ,C (39) 

EQUIVALENCE ( U (1) ,WE ( 1) ) , (V ( 1) f W N ( 1) ) 

READ (5,40) 18 ,JM,IMH, JHH 

40 FORMAT (4 15) 


IMP1=IM+1 


JMP1=Jtt+ 1 
C CALL SHAPE1 

CALL SHAPE 1 
DO 548 J=1,MM 
DO 549 1= 1 , I MP1 
F (I , J ) - 0 . 0 
G (I , J) -0 . 0 
U (t , J) = 0 . 0 
V (I , J) -0 . 0 
XU (I, J) =0.0 
i, XV(I,J)=0.0 

, 549 CONTINUE 

V 548 CONTINUE 
1 REA D (5, IS) Y , DX, DY 

18 FORMAT (3 FI 0.0) 

P=3. 1415926/180.0 

DX= DX *P 

MK= JM+2 

DO 100 1=1, mm 

AR= ( Y+D Y*FLOA T (I) ) *P 

F1 = CCS(AR) 

X(I)=DX*F1 
C (I ) =F 1 

BET (T)=F1*X (I)*14.58E“5 
100 CONTINUE 
Y=D Y* P 

CD=1 4 .58/ (DX* Y) 

DO 9997 J= 1 , J MP 1 
XM(J) = (X (J) + X (J + 1) ) /2.0 
FJU ( J) =CD* (X(J) -X (J + 1) ) 

9997 CONTINUE 

READ (5, 16) A H , H , REL , DT , MAT 

16 FORMAT (4 El 5. 7 ,15) 

WRITE (6, 17) AH, H, REL 

17 FORM AT ( IH , 1 5 HEDD Y VISS OSIT Y=E1 5. 7, 1 0 X,6HDEPTH=.E 1 5. 7 ,1 0 X, 
1 1 5H RELAX ACCURACY=E15.7) 

• I READ (5,10) IT 

10 FORM AT (1 5) 

GOTO (1,2), IT 
1 CONTINUE 
C CALL WIND 

CALL WIND 
P=0. 01/4. 0E5 
DO 540 J = JMH , JMP1 
y MJ - X M (J) 


XM1=XM (JM1) 

DO 54 1 I=IMH, IM 



TW(I,J) = (XMJ* (-ME (T,J)-WE (IM1, J) } +XM1*(WE (I # JH1) +WE(IM1,JM1) ) 
1 + Y*(WN{I,J)+WN{I,JM1)-WN(IF J 1,J)-WN{II']1,JM1)))*P 

1 CONTINUE 

T 5? ( 1 , J) =TW (IMM1, J) 


TW (2,0) =TW (IM r J) 

TM(IMP1 , J) = TW (IHH, J) 

540 CONTINUE 
GO TO 3 

2 CONTINUE 

C CALI CURL 

CALL CURL 

READ (5, 60 04) (TW (1,0) ,I = 3,IM) 
6004 FORMAT (5E15, 7) 

TW (IMP1 , J) =TW (3, J) 

TW (2 , J) =TW (IM , J) 

TW( 1, J) =TW(IWM1,J) 

3 CONTINUE 

C INITIAL GUESS 

DO 4000 J=4,JM 
CJ=1 . 0/ (C (J) *Y* 14. 58E- 5} /4 . 0 


L J=L (J) 

DO 3 6 0 LL- 1 , LU 
11= K A (LL, J) 

IG= KB (LL, J) 


MP= II +1 G 

362 DO 361 1= II , I G 
m~n P-I 
K=fi M + 1 


P=P (K f J) - (TW(MM, J) +TW (K,J) ) *CJ 
? {MM, J) -P 
361 CONTINUE 

IF (II - NE . 3) GO TO 360 
F(IMP1, J) =F (3 , J ) 

360 CONTINUE 

F (1 , J) =F(IMH1 ,J) 

F (2, J) = F (IM, J) 

4000 CONTINUE 
C INITIAL GUESS 


CCC=2.0*Y 
DO 4049 0=3, JM 
JP1 = J+1 


LJ-LV (J) 

DO 405B L1=1,LJ 
11= LA (LL, J) 

IG— LB (LL ,3) 

DO 4053 I=II , IG 


IP1=I*1 


p = (F (I, J) +F (IP1,J)-F (I, JP1)-F (I PI ,JP1 ) ) /CCC 
0(1. J)=P 
XU (I , J) =P 


P= (F {IP1 ,3) +F (I PI , JP1 ) “F (I, J) -F (I,JP1 )) /XMJ 
V(I.J)=P 




XV (I , J) = P 
3 CONTINUE 
'58 CONTINUE 


XU (IMP1 , J) =XU (3 , J) 
XV(IMP1 , J ) = XV (3 , J ) 
XU ( 2 , J) =X U (IM ,3) 
XV(2,J)=XV(I«,J) 

XU (1 ,0) = XU (IM M 1 , J } 
XV(1,J) =XV(IMM1,J) 
T! (I Fi P 1 , J) =0(3, J) 



V (IMPl, J) =V (.-5,0) . 

u {2 # j)=u ( in r j) 

' H V (2, J)=V (IM,J) 

" ► ' o (1 „ J)=U (IMMl ,3) 

V (1 ,3) = V (IMMl , J) 

ICQ 49 CONTINUE 

AR=AH/6. 37EB/6-37E8 

ARY = AR/Y 

CC Y = 2 . 0 * Y 

Y4 = Y *0. 25 

RETURN 

END 

SUBROUTINE SHAPE1 

C SUBROUTINE TO BE USED FOR THE BAROTROPIC PART OP THE WORLD OCEAN 

C CIRCULATION WHEN THE SHAPE OF THE OCEAN IS CHANGED. 

C DATA; LY, LA AND LB FOR 15 Q Q. 

C KM, M1(K),M2(K) AND M3 (K) FOR MPP. 

C LV r LA AND LB ARE THE SAME AS THOSE USED FOR THE SUBROUTINE 

C SHAPE 2 (BAROCLINIC PART) , BUT THE DATA FOR MPP ARE DIFFERENT 

C FROM THOSE USED FOR THE SHAPE2. 

C LX1, LX2 , IY1, LY2 WHICH SPECIFY THE AREAS INCLUDING THE SEPARATE 

C ISLANDS. 

COMMON AH,DT,Y,H,RE1,XM(3 9),FJ0 (39), X (3 9) , F (147,39) , E (147,39) , 

10(1 47,3 9) ,V (147,3 9) ,G (147 , 39) ,X(J( 147, 39) , XVC147, 39) ,TH< 147, 39) , CD , 

23 ET (39) ,Y4,CCY,AR,ARY,RES (3) r DDT, 

3 MAT, I (3 9) , L V (39) ,LA (10,37) ,LB (1 0, 37) , KA (1 0,37) , KB (1 0 ,37) , MATS , 

4MPP (147,3 9) ,MQQ (147 ,3 9) , INDIN (3) , IND{3) , NRES (3) , 

5LX1 (3) , LX2 (3) ,1Y1 (3) , LY2 (3) ,IM, JM,IHP 1, JMP1, IMM 1 
IMM 1-IM- 1 
I MP 1 =IM+ 1 

JMP1 = JM+ 1 
J MH= 3 
I MH =3 

DO 30 3 = 1 , 3 M 
DO 31 1= 1 , I MP 1 
MQQ 

31 CONTINUE 
30 CONTINUE 

DO 32 J = 3 , J M M 1 
READ (5,3 3) M 

33 FORMAT (12) 

READ (5, 3 4) (LA (I, J) , LB (I,J) ,1 = 1, M) 

34 FORMAT (7 (214, 2X) ) ; 

DO 35 1=1 ,M 

11= LA (I, J) 

I G=L B (1,0) 

DO 36 N = II , IG 
MQQ (N,J) =1 
36 CONTINUE 

35 CONTINUE 
LV (J) =M 

HQQ(1,0)=MQQ(I»M1,J) 

MQQ(2,J)=MQQ(IM,J) 

MQQ(IMP1 ,3)=HQQ (3 , J) 

32 CONTINUE 

DO 3006 3=3,38 

« READ (5, 3007) MM 
7 FOR MAT (I 2) 

READ (5, 3 0 08) (Ml ( K) ,M2 (K) , 153 (K) # K=1, MM) 

FORMAT (8 (214,12)) 

DO 3009 K=1,MM 
K5~ M 1 (K) 

K6- M2 (K) 

K7= M 3 (K) 

DO 3010 I=K5, Kf> 


3008 



MPP 1 {I, J) =K7 
v BO 1 0 CONTINUE 
r 3G09 CONTINUE 
'300 6 CONTINUE 

DO 3011 J ~ 3 , 3 8 
READ (5, 3007) MM 
LI (J) =MM 

READ (5,3012) { KA 1 (LL,J) ,KB1 (LL, J) , 1I,= 1 , MM) 

3 012 FORMAT (1615) 

3011 CONTINUE 
DO 6 003 

MPP 1 (1, J) =MPP 1 (IHM1,J) 

MPP1 (2, J) =MPP 1 (in, J) 

MPP1 (IMPl,J) = HPP1 (3,J) 

6003 CONTINUE 

DO 6 0 ON J = 3 , J M 

JP1=J+1 

UM 1= J-1 

DO 6005 1=3,1 M 

IF (MPP1 (I, J) . GT. 1 ) GO TO 6005 - 
IF (MPP1 (1-1, J) -EQ..2) GO TO 6005 
IF (HPP1 (I, J?1) -EQ.2) GO TO 6005 
IF (KPP1 (I, Jill ) . EQ, 2) GO TO 6005 
IF (M.PP1 ( 1 + 1 , J) . EQ.2) GO TO 6 0 05 
KPP1 (1,3 ) =0 
6005 CONTINUE 

HPP1 (IMP1,J)=MPP1 (3,J) 

MPP 1 (2, J) = MPP 1 (TM , J) 

MPP1 {1, J)=MPP1 (IHM1 ,J) 

6004 CONTINUE 

WRITE (6 , 6016) (J,J=3,38) 

6016 FORM AT ( 1 H ,5X,3 8I3/) 

DO 6012 1= 1 , I MP 1 

WRITE (6, 6014) I, (MPPl (I , J) ,0=3,38) 

WRITE {6, 6015) (HQ Q (I , J) , J- 3, 3 8) 

6012 CONTINUE 

6014 FORM AT ( 1 H ,13, 2X, 3813) 

6015 FORMAT ( 1 H ,6X,38I3) 

READ (5,2) (IX 1 (K) ,LX2 (K) , LY 1 (K) ,LY2 (K) ,K=1,3) 

2 FORMAT (4 15) 

RETURN 

END 

SUBROUTINE WIND 

COMMON AH,DT, Y, H, REL,XM (3 9) , FJU (39) ,X (39) ,F (147,39) ,E(147,39) , 

1U (1 47,3 9) , V (1 47,3 9) ,G (1 47 ,39) ,XU ( 147, 39) ,XV (1 47, 39} ,TW ( 147, 39) ,CD 

2 BET (39) , Y4,CCY,AR , ARY, RES (3) , DDT, 

3 HAT, 1(39) ,LV (39) ,LA (1 0, 37) ,LB (1 0, 37) , KA (10,37) ,KB (10,37) ,MATS, 

4 MPP (147,39) , MQQ (1 47,3 9) , INDIN (3) ,IND (3) , NRES (3) , 

5LX1 (3) ,LX2 (3) ,LY1 (3) , LY 2 (3 ) , I H, JM , I MP 1, 3 MP 1 , IMM1 

DIMENSION WE( 147,7) ,WN ( 147,7) 

RETURN 

END 

SUBROUTINE CURL 

RETURN 

END 






EAROTROPIC CIRCULATION IN A MULT I- CONNECTED GLOBAL OCEAN 

MPP=0 IF TEE FOUR NEAR ESR SURROUNDING POINTS ARE WATER POINTS, 

=1 IF ONE OF THE FOUR NEAREST SURROUNDING POINTS IS ON THE 
COAST OF THE MAIN LANDMASS (AFRICA, EUROPE, ASIA), 

=3 INSIDE THE ANTARCTICA, 

=30 ON THE COAST OF THE ANTARCTICA, 

=4 INSIDE AUSTRALIA, 

=40 ON THE COAST OF AUSTRALIA, 

= 5 INSIDE NEW ZEALAND, 

=50 ON THE COAST OF NEW ZEALAND. 

MULTIPLY THE STREAM FUNCTION VALUES BY R*H (B^EARTH* S 
RADIUS, H=OCEA N DEPTH) TO GET THE STREAM FUNCTION VALUES 
FOR THE VOLUME TRANSPORT (INTEGRAL OF THE HORIZONTAL VELOCITY 
OVER THE WHOLE DEPTH. 

MQQ=1 FOR THE WATER POINTS, 

=2 FOR THE MAIN LAND MASS, 

=3 FOR THE ANTARCTICA, 

=4 FOR AUSTRALIA, 

=5 FOR NEWZEALAND. 

TEE SUBROUTINE HAJIME SPECIFY THE SHAPE OF THE OCEAN (MPP AND 
MQQ) , THE WIND STRESS TO COMPUTE ITS CURL AND OTHER CONSTANTS 
AS WELL AS THE INITIAL FIELD OF THE STREAM FUNCTION. 

THE SUBROUTINE .RES ID GIVES THE CHANGE OF THE SUM OF THE 
RESIDUALS (RES) AT EACH SEPARATE LAND MASS WHICH TAKES PLACE 
BY THE UNIT INCREASE OF THE STREAM FUNCTION AT EACH SEAPARATE 
LAND MASS POINT, JUST AS IN CASE OF THE BLOCK RELAXATION. 

F AND G DENOTE THE STREAM FUNCTION AT TIME STEP N+1, N, 
RESPECTIVELY, XU,XV,U, A ND V THE BAROTROPIC COMPONENTS AT TIME 
STEPS N+1, N, RESPECTIVELY. 

WIND STRESS CURL=- TW/ ( 2 . 0* Y* X ( J) ) /R*H 
X , XM , Y =GRID SIZE 

A H=COEFFIC IENT OF EDDY VISCOSITY 
DT IS THE TIME STEP. 

H IS THE OCEAN DEPTH- CD IS THE TIME (DAY). 

E IS THE WORKING SPACE FOR RELAXATION. 

IM DENOTES THE EASTERNMOST POINTS TO BE SWEPT. (EX- 146) 

JM DENOTES THE NORTHERNMOST PONTS TO BE SWEPT. 

THE BACKWARD DIFFERENCING SCHEME IN TIME APPLIES EVERY MAT 
TIME STEPS. MATS CONTROLS THE SWITCHING FROM THE LEAP-FROG TO 
BACKWARD OF FROM THE BACKWARD TO LEAP-FROG SCHEMES. 

PRINTING IS DONE EVERY NPRINT TIME STEPS. 

WRINTING IN MAGNETIC TAPE IS DONE EVERY NTAPE TIME STEPS. 

WHEN (IN) EQUALS ONE, THE IMAG-TH RECORD IN THE MAGNETIC TAPE 
IS PICKED UP, AND THE INTEGRATION GOES ON UNTIL NN=NEND. 

REL= RELAX AT ION ACCURACY 

WIND STRESS CURL=-TW (I,J)*H (DEPTH }/{2.0*Y*X (J) *RAD) RAD= 
RADIUS OF THE EARTH) 

I ND (1) =30 ANT ARCTIC COAST 
IND (2 ) =40 AUSTRALIA COAST 
IND { 3) =50 NEWZEALAND COAST 
INDIN (1) =3 ANTARCTIC INLAND 
INDIN (2) =4 AUSTRALIA INLAND . 

INDIN (3) =5 NEWZEALAND INLAND 

FIGURE 1 EXPLAINS THE MEANING OF LX1,LX2, IY 1 AND LY 2. 
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FIG. 1 

L,KA,KB SPECIFY THE STREAK FUNCTION POINTS TO BE SWEPT AND 
OTHER LAND POINTS TO BE SKIPPED. 

FIGURE 2 EXPLAINS THEIR MEANINGS. 

* POINT TO BE SWEPT 


1 = 4 KA (1,J)***KB(1,J) KA(2,J) ** ***K B { 2, 3) KA (3, J) ****KB{3,J) 

L (J) =3 


C ' FIG. 2 

C 

c 

C THE COMPUTATION STOPS IF THE NFRGY EXCEEDS A CERTAIN VALUE, OR 

C IF THE NUMBER OF ITERATION IN RELAXATION EXCEEDS A CERTAIN VALUE. 

C THESE VALUES ARE SPECIFIED IN THE PROGRAM. THEY SHOULD BE CHANGED 

C IF NECESSARY. 

C THE COEFFICIENT OF THE OVER-RELAXATION IS GIVEN AT THE BEGIN- 

C NING. 

C 

C 

COMMON AH,DT, Y,H,REL,XK (39) , FJU (39) ,X (39) ,F (147,39) , E (147,39) , 

10 (1 47,39) ,V (1 47,3 9) ,G (147,39) ,XU{ 147, 39) , XV < 1 47, 3 9) ,T W < 147, 39) , CD, 

2 BET < 3 9) , Y4,CCY, AR , ARY, RES (3) , DDT, 

3ttAT,L(3 9) ,LV (39) , LA (10, 37) , LB (10,37) , KA (1 0 , 37 ) , KB ( 1 0 ,37 ) ,BATS, 

4 MPP (147,3 9) ,MQQ (1 47,3 9) , INDIN (3) , IND (3) , NRES (3) , 

51X1 (3) ,LX2 (3) ,LY1 (3) , LY 2 (3) , IH, JM,IHP 1, JMP 1, IH8 1 
" SOR = 1 .77 5 
SOR = 1 *87 5 
SOR =1.825 
SOR= 1.85 
SOR 1 = SOR — 1.0 
READ (5,10 00) IN 

1000 FORMAT (15) 

IF (IN.EQ. 1) GO TO '1001 

CALL HAJIME 

CALL RES ID 

CD= 0.0 

NN=0 

MAG = 0 

GO TO 5555 

1001 CONTINUE 

READ (5, 1000) IMAG 
IF (I MAG . EQ. 1) GO TO 77 
DO 7 J=2 , I MAG 
READ (10) 

7 CONTINUE 

# 7 READ (10) XU, XV,F,G,U,V,TH,MPP,HQQ,L, KA, K B, LV ,LA,LB, AH, DT, MAT, 

1 Y,REI,XM,FJU, X,BET,LX1,1X2,LY 1, LY2, IND, INDIN, NH, REND, RES, 

2NRES, ARY, AR,Y4,E,CCY, MATS ,CD ,IH , JM, IM PI , J MP1 ,IHM 1 , MAG 
5555 CONTINUE 
A H - 1 . OE 9 
AH= 3 . 0E8 

AR- AH/6. 37E8/6. 37 E8 
ARY=AR/Y 
I HAG=NN+ 1 


&SN-U,U 


READ (5,5) REND, HAT, HP HI NT , NT APE, DT, PEL 
*, 5 FORMAT (4 15,2515.7) 

" K RRITE(6,6) IMAG,flEND, DT,REL,AH, MAT, NT APE, SPRINT 

^ 6 FORM AT ( 1 H , 16HTIME STEP FROM I5,2X,3HT0 15 , 2 X, 3HBT=E 15 „ 7 , 2X, 4H REL 

1=E15*7,2X,3HAH=E15-7,2X,4HMAT=I3/1H , 5HT APE=I4, 5X , 6 HPRIRT=I 4) 


£0 DAY=DT/0, 864E5 
MQQ (4 9, 5) =1 
MQQ (50,5) =1 
MQQ (49, 3) =3 
MQQ (50, 3) =3 
39 NN = NN + 1 


CD=CD+DAY 

N NN= (NN/M AT) *11 AT 

IF( NNN.EQ, NN) GO TO 25 

MAT S = 0 

DDT=DT 

DO 1 J= 3 , JM 

DO 2 1=1, IMP 1 

R=G (I ,J) 

G(I,0)=F (I,J) 

F(I, J)=R 
2 CONTINUE 
1 CONTINUE 
GO TO 26 

25 WATS=1 
DDT= DT* 0 . 5 
DO 42 0=3 , JM 
DO 43 1=1 ,IMP1 
U (I, 0) =XU (I, J) 

V (I, J)=XV (1,0) 

G(I, J)=F(I,J) 

43 CONTINUE 

42 CONTINUE 

26 CONTINUE 

DO 30 0=3, JM 
0P1=J+1 

OM 1=0-1 

0 M2 =0-2 
XMJ = XM(J) 

XM1 = XM ( J Ml ) 

C 1= Y/XMO 

S1-Y/XM1 

CB= X (OP 1) /Y 

CS=X (0) / Y 

B=X ( JM1 ) /Y 

BET J=EET (0) *CCY 

LO=L (0) 

DO 31 LL=1,LJ 
11= K A (LL, J) 

MM=KB (LL, J) 

DO 310 I = II,MM 
I P 1 = I + 1 
IHl=I~1 
IM 2= 1-2 
UIJ=U (1,0) 

VIO = V (1,0) 

UP10 = U(IP1,J) 

DMl 0=U ( I M 1 , J) 

UIP 1 =U (I , 0P1 ) 

UM1P1=0 (IM1, JP1) 

DIM 1=0(1 ,0M1) 

UMl M 1 =U (IM1/JM1) 

UP1M 1=0 ( IP 1 , JM 1) 

UIM2 = U(I, JM2) 

VT 0 2 = V (I , 0M2) 


UK2J=U (1112, J) 

VM2 J- V (I M2 , J) 

UM1M2-U (IM1 , JM2) 

, VM1M2=V (IM1,JM2) 

U M2 K 1 =U { I M2 , J M 1 ) 

" VM2H1=V (IM2,JM1) 

VP1J=V(IP1,J) 

VM1JJ = V(IM1,J) 

VIP 1=V (I,,1P1) 

VIM 1=V(I, Jfil) 

VK1P1=V (I Ml , JP1) 

VM1 M 1=V (IK1/JM1) 

VP 1 M 1=V (IP 1 , J M 1) 

MP= MFF (I , J) 

IF<MP.EQ. 0) GO TO 304 

IF (MP.EQ. 30) GO TO 304 

IF ( MP-EQ. SO) GO TO 304 

IF(MQQ (IP1, J) . NE. 2) GO TO 27 

UP1J=-UIJ 

VP1J=-VIJ 

27 IF(MQQ(I, JP1) .NE. 2) GO TO 271 
TJIP 1 =-UI J 

VIPl=~VIJ 

271 IF (MQQ<IM1, OP 1) . NE.2) GO TO 272 
U M 1 F 1 -- D M 1 «1 

VM1 P1*-VH1J 

272 IF (MQQ (IM2 # J) .NE. 2) GO TO 28 
0 M2 J=— U Ml J 

VM2J*-VM1 J 

28 IF ( MQQ ( I M 2 , JM 1 ) • NE. 2) GO TO 29 
U M2 M 1 =- U M 1 M 1 

VM2M1=-VM1M1 

29 IF(MQQ{IM1,JM2) . NE.2) GO TO 291 
UM1 M2--UM1M1 

V M 1 M 2»“ V M 1 Ml 

291 IF(MQQ (I f JM2)-. NE. 2) GO TO 292 
UIM2=-UIM1 

VIM2«-VIM1 

292 IF {MQQ(IP1, JMl) . NE.2) GO TO 303 
0P1 Ml ——til Ml 

VP1M1=-VIM1 

303 CONTINUE 

304 DF--BET J * (XV (I f J)+XV(IMl,J)+XV (IM1 r 0Ml) + XV (I, JMl)) *0.25 

S S- (C 1*{ UM 10 + UI J-UM2J-UP1 J) - (UTP1-UIJ + UM1 P1-UM1 J) *CB42. 0*CS* 

1 (DIJ-UIMUUM1 J-UMTH1) 4S1* (UP1M1 +UM2K1 -UIM1-UM 1M1 ) - (UIM1-UIM2+UM 1M1 
2-UMlM2)*3 ) *ARY + TV! (I, J) +DF 
DF=CS* (VIJ-VI M1-VM1J4VK1M 1) 

SS= AR/XM J* (C 1* (-3. 0* (VI J-VM1 J) + VP1J-VM2J)-DF 

14 (VIP1- VIJ-VM 1P1 + VM1J) *CB) 4AE/XM1 * (SI * (3. 0* (VM1H1-VI81) +VP1M1 
2-VM2M1) 4- (VM1M1-VM1M2-VIM1 +VIM2) *B +DF) +SS 
4 E (I, J) =- SS*DDT 
310 CONTINUE 
31 CONTINUE 


E (2,a) = E (IP5 , J ) 

E (IMP 1,0) =E (3, J) 
30 CONTINUE 
IP=0 



7 IP=IP4l 
** ISt, ANUS 


M AR = 0 


c *****antapctica 

C ***** N £ W ze ALAND 
c * ♦ ***& USTR ALIA 

DO 5021 K«1,3 
RESK-RES (K) 

M P" I I'f D (K) 



1J F i J F — X iM JJ 1 \ l\ } 

11= LX 1 (K) 

IG=LX2(K) 

‘ h I M2 =L Y1 (K) 

* " JM2=LY2(K) 

R = 0. 0 

DO 5016 J=IM2,JM2 

DO 5015 I=II,IG 

IF(MPP(I, J) -NE- HP) GO TO 5015 

R=R * E (I , J) 

5015 CONTINUE ' 

5016 CONTINUE 

IE ( A BS ( R J.LE.REL) GO TO 5021 

DF=~ R/R ESK 

MAR=MAR-M 

IF(K-EQ.I) IM 2=2 

DO 5022 J=IM2 # JM2 

JP1=J+1 

J N 1 = J - 1 

CS= Y/ X { J) 

CB= 2-0* (CS+1.0/CS) 

SS=X (JP 1 ) /Y 
C1 = X (JM1 ) /Y 
DO 5023 I=II,IG 
HPPI J=KPP (I, J) 

HMMM= (MPPIJ-MP) * (MPPIJ-HPMP) 

IF(MMMH.NE.O) GO TO 5023 

F (I , J) =F (1,0) + DF 

IM1 = I-1 

IP1=I+1 

CCC=DF*CS 

E (I, JP1 ) = E (I, JP1) +SS*DF 
E {I, J«1) =E (I r JK1) + C 1* DF 
E (IP1, J) = E (TP1, J) +CCC 
E (I M 1 , J) =E (I Ml, J) +CCC 
E (I, J) = E (I, J) -CB* DF 

5023 CONTINUE 
5022 CONTINUE 

IF(K-EQ-2) GO TO 5021 
IF (K . FQ- 1) GO TO 5024 
DO 5 026 J = 8 , 1 2 
E (2, J)=E (IN, J) 

.5026 CONTINUE 

GO TO 5021 

5024 CONTINUE 

E (in PI, 5) =E (3,5) . . 

E(2,5)=E(IH,5) 

CCC= DF* Y/X (4) 

E (3,4) =E (3,4) +CCC 
E (I 4) = E (IM, 4) + CCC 
E(in?1,4)=E (3,4) 

E (2 , 4) = E (IM, 4) 

5021 CONTINUE 
c *****newzeai.and 

C *****ANT ARCTICA 

DO 32 J=4,JM 
J«1=J-1 
JP1=J+1 
X0=X (J) 

CS= Y/XJ 
SS=X (0P1)/Y 

ci=x (in i) /i 

CB-2 . 0* (CS+ XJ/Y) 

LO=L (0) 

DO 33 LL = 1 ,L J 
11= KA (LL, J) 



a 



DO 330 I=II,MM 

, IF(MPP(I,J) .GE.3) GO TO 330 
► K=E(I,J) 

IF (ABS(R) .LE. REL) GO TO 330 

MAR=MAR+1 

DF=R/CB*SOR 

E{I, J)=-R*SOR1 

F (I, J)=F (1,3) + DF 

IP1=I+1 

CCC = D P*C S 

E (IPl ,3) = E (IP 1 ,3) +CCC 
IM 1=1-1 

E (IM1,J) =F, (IH 1,3) +CCC 
E (I, JP1) = E (I,JP1) +SS*DF 
E (I , JM1) =E (I, JM1) +C1*DF 
IF(I.NE.3) GO TO 331 
E{IMP1, JP1)~E (3, JP 1) 

E (IMP 1, JH1) =E (3, JM1) 

E (I MP1 ,3} = E (3 , J) 

E (IM, J) =E <2, J) 

GO TO 330 

331 IF(T.NE.IM) GO TO 332 
E (2, J)=E (IM, J) 

E (2, 3 PI) - E ( IM , JP 1 ) 

E (2, JM1) = E (IM , J M 1 ) 

E (3 ,J)=E (IMP1 ,3) 

GO TO 330 

332 IF(I.NE.4) GO TO 333 
E (IMP1 ,3) =E (3,3) 

GO TO 330 

333 IF(I.NE.IMMI) GO TO 330 
E (2 , J) -E (IK, 3) 

330 CONTINUE 
33 CONTINUE 
32 CONTINUE 

IF(MAR.EQ.O) GO TO 36 
" IF (IP. GT. 1000) GO TO 1112 
IF (IP.NE. 1) GO TO 37 
WRITE (6,38) IP, MAR 

38 FORMAT ( 1 H , 7X , 6HKA ISU=T 5, 5X, 12HTEN NO. KA2U=I5) 

GO TO 37 
- 36 EN= 0 . 0 

WRITE (6, 38) IP, MAR 
DO 7 0 3=3,311 
F (1/ J)=F (IMM1 ,3) 

F (2,3) =F ( IM , 3 ) 

F (IMP 1,3) =F (3,3) 

70 CONTINUE 

DO 23 3=3, 3M 
JM1=J~1 
JP 1=3 + 1 
X MJ = X K ( J) *2.0 
DF= XM J* Y4 
LV3 = L V (3) 

DO 24 13=1, LV3 
11= LA (13,3) 

IG=LB (LJ, 3) 

DO 244 1= It , I G 
IP1=I+1 

U (I , 3) = XU (I, J) 

V (I,J)=XV (1,3) 

XU (1,3) = 

1 (F (I, J) +F (IP 1,3} -F (I, JP1 )-F (IP 1, JP1) ) /CCY 

XV (I ,J) = 

1 (F (IP1 ,3) +F {IP1, JP1) -F (I, J)-F {I, JP1) ) /XMJ 



x>N- juH (AlS 11,0} '^2 +aV 11 ,0 j « •/ J 

244 CONTINUE 
«. 24 CONTINUE 
„ XU (IMP1, J) =XU (3 r J) 

XV (IMP1,J) -XV ( 3 , J ) 

U(IMP1,J)=U(3,J) 
y (IMP 1, J) = V (3 , J) 

XU (2,0) =XU (IN, J) 

XV{2, J) = XV (IK ,J) 

XU (1 ,0) = XU (IMM1 , J) 

XV ( 1, J) -XV (IN M 1 , J ) 

U (2, J) =U ( I M , J ) 

V (2 , J) = V (IM,J) 

U (1 , J)=U (IMH1 r J) 

V (1 , J)=V (IMM1 , J) 

23 CONTINUE 

WRITE (6, 34) NN,CD,EN,F (1,3) 

3 4 FORM AT (If! , 5H ST EP -15, 5X ,5 HTI ME=E 1 5. 7 , 5X,7HENERGY=E15.7,5X,4HAC 
1 ,4) 

IF(EN.GT. 5.0E1) GO TO 1111 

IF (NN.GE. NEND) GO TO 88 

IF { MATS* NE.O) GO TO 10 

IF(NN.NE. ( (NM/NTAPE) # NT AP E) ) GO TO 9 

MAG-N AG+ 1 

WRITE (10) XU, XV # F,G,U,V,r W, MPP, MQQ, L, KA, KB, L V ,L A, LB, AH , DT, M AT, 
1Y,R2L,XM,FJU, X,BET,LX1, LX 2 ,LY1,LY 2,IND,INDIN, NN , NEND, RES, 
2NRES, ARY, AR, Y4,E,CCY, MATS , CD ,IM , J N, IM PI , J MP 1 , 1 SM 1 , NAG 
WRITE (6,8) MAG 

9 IF ( NN , EQ. ( (NN/NPRI NT) * N PR I NT) ) GO TO 88 
10 ENN = EN 

IF ( H ATS. EQ. 0) GO TO 39 
MAT S = MAT S + 1 

IF (MATS. EQ. 3) GO TO 4 6 
DO 4 4 J = 3 , J M 
DO 45 1=1 ,IMP1 
F (I, J) =G (I,J) 

U (I, J) = XU (I, J) 

" V(I,J)=XV(T,J) 

45 CONTINUE 
44 CONTINUE 

GO TO 26 

46 DO 47 J=3,JM 
DO 48 1=1 ,IMP1 
R = F (I, J) 

F(I, J)=G <I,J) 

G (I , U ) = R 
48 CONTINUE 

47 CONTINUE 

DO 49 J = 3 , «7 M 

JP1=J+1 

XMJ = XM (J) *2.0 

LVJ = LV(J) 

DO 55 L J= 1 ,LV J 
II=IA (I J , J) 

IG=1B(LJ,J) 

DO 550 I = II ,1 G , * 

IPl-I+1 

0(I#J) = 

1 (F (I,J) +F{IP1, J) -F {I,UP1) -F{IP 1, JP 1) )/CCY 

V (I, J) ” 

1 (F(IP1,J)+F(IP1,JP1)~F (I,J)-F(I,JP1) )/XMJ 

550 CONTINUE 
55 CONTINUE 

U (1 , J) = U (IMM1 ,0) 

V (1 , J) = V (IMH1 , J) 

U (2, J) =U (IM, J) 


V yl r 0) -V l.Lw,*J) . 

U (IHP1,J)=U (3 ,0) 

V (IKP1, J) =V(3,J) ; 

*' ^49 CONTINUE 

* MATS=0 
NN= NN+1 
DDT-DT 
NNT= NNT + 1 
CE=CD+DAY 
GO TO 39 
88 K = 1 
MP= - 1 

613 MP=KP+1 
MM=K~MP 
IG= 1 8*MP + 3 

IF(MM.EQ.O) GO TO 614 

II=TG+17 

GO TO 615 

614 11=3 8 

615 WRITE (6, 7 01} 

701 FORM AT ( 1H0, 15HSTR EAM FUNCTION) 

DO 7 02 J= 1 , IM Pi 

WRITE (6,7 0 3) J, (F (0,1) , I=IG,II) 

703 FOR M AT ( 1 R , 1 3 , IX, 18F7 - 4 ) ; 

702 CONTINUE > 

IF(MM.NE.O) GO TO 613 

IF (NN.LT.NEND) GO TO 10 
M AG=MAG+ 1 

WRITE (10) XU, XV,F ,G , U, V ,T W , KPP, KQQ, L, KA , KB,LV ,LA,1B, AH, DT, MAT, 

1 Y ,R EL,XM , FJU, X, BET , LX 1 , LX 2 ,IY 1 , LY 2 , IN D, INDI N , NN , NEND ,R ES , 

2NRES r ARY, AR,Y4 ,E,CCY, MATS , CD , I M , J M, IM PI , JMP1 , IMM1 , M AG 
WRITE (6,8) MAG 

8 FORK AT( 1H0,32HNO. .OF RECORDS IN MAGNETIC TAPE=I3) 

GO TO 1111 

1112 WRITE (6,38) IP,MAR 
1111 STOP 
END 

"" SUBROUTINE RESID 
RETURN 
END 

SUBROUTINE HA 01 ME 

RETURN 

END 

SUBROUTINE RESID 

COMMON AH,DT,Y,H, REL, X K (3 9) , FJU (3 9) ,X (39) , F (1 47 ,3 9) ,E ( 1 47, 39) , 

1U (1 4 7,39) ,V (1 47,3 9) ,G (147 ,39) ,XU( 147, 39) , XV ( 1 47, 3 9) ,TW ( 147, 39) ,CD, 
2BET (39) , Y4,CCY,AR,ARY,RES (3) , DDT, 

3 MAT, I (3 9) ,LV (39) ,LA (10, 37) , LB (1 0, 37) , KA (10,37) , KE (10,37) , MATS, 

4MPP (1 47, 3 9) ,MQQ(1 47,39) , INDIN (3) , IND (3) , NRES (3) , 

5LX1 (3),LX2 (3) ,LY1 (3),LY2(3) , IM, JM,IMP 1 , JMP1 ,IMM1 
DO 1 K= 1 , 3 
DO 8 J=1 , JMP1 
DO 9 1= 1 , IMP 1 
E (I , 0) =0. 0 

9 CONTINUE 
8 CONTINUE 

S=0. 0 1 

I 1=LX 1 (K) 

I2=LX2(K) 

J 1 = L Y 1 ( K) 

J 2=L Y2 (K) 

MIN=IND(K) , 

MIN IN = INDIN (K) 

NUM = 0 

IF(K.EQ.I) 01=2 
DO 3 0=01,02 



non 


v)M5-d-J 
JP1=.1 + 1 
CS= Y/X ( J) 

„ C1 = X (0M1)/Y 
^ SS=X(UP1)/Y 

C8=-2«0* (CS+1 .0/CS) 

DO 4 1=11,12 
MP=MPP(I,J) 

MMM= (MP-MIN) * (HE- MI NIK) 

IF (MMM* NE. 0) GO TO 4 
IP 1 = 1 + 1 
IM 1=1-1 
NUM= NUM + 1 
E (I, J)=E (I, J) +CB 
E(IP1,J) = E (IP1, J) +CS 
E (IM 1, J) ~E (IM 1, J) + CS 
E (I, JP1) = E(I, JP1) +SS 
E (I, JM1) =E (I , JM1) +C1 

4 CONTINUE 
3 CONTINUE 

IF (K . GE. 2) GO TO 10 
J 1= 3 

DO 11 J = 3 , 3 
CS= Y/X (0) 

E (3, J) = E (3, J) +CS 
E (IM, J) =E (IM, J) +CS 
11 CONTINUE 
10 CONTINUE 

DO 5 0= J 1 , 02 
DO 6 1=11,12 

IF (MPF(I, J) .EQ. MIN) S=S+E(I,J) 

6 CONTINUE 

5 CONTINUE 
NBES (K) = NUM 
RES (K)=S 

WRITE (6, 7) K , NUM, S 

7 FOR MAT ( 1 f! ,2 H K=I2 , 2 X, 1 2 HG RID POINTS=I 3, 2X ,4 HRFS = E1 5 - 7) 
1 CONTINUE 

DO 16 J=1,18 
DO 17 1=1, IMP 1 
IF (MQQ (I, J) -EQ. 1) GO TO 17 
• MP=MPP (I , J) 

IF (MP.EQ. 2) GO TO 17 
IF(MP.EQ, ( (MP/3)*3) ) GO TO 18 
JF(MP.EQ. ( (MP/4) *4) ) GO TO 19 
MQQ (I,J) =5 
GO TO 17 

18 MQQ (I,J) =3 
GO TO 17 

19 HQQ (I f J)=4 
17 CONTINUE 
16 CONTINUE 

WRITE (6, 15) (1,1=1, 38) 

15 FORM AT (1H ,4X,38I3) 

DO 12 1=1 ,IMP1 

WRIT E (6 , 14) I, (MPP(I, J) ,J = 1,33) 

WRITE (6,1 3) I, (MQQ{I,J) ,J=1,38) 

• 12 CONTINUE 

13 FORMAT { 1H ,I3,3X,38I3) 

14 FOR M AT ( 1 H ,I3,2X,38I3) 

RETURN 

END 

SUBROUTINE HAJIME 


4S^ 


DATA: 

IM, JMH 



U UCJSAN StfAPb (bUbHUUXi.NL bHfiPi. J) 

C X , DX ,DY (DX AND DY ARE TEE GRID SIZES IN LONGITUDE AND 

LATITUDE (DEGREES), RESPECTIVELY. Y IS THE LATITUDE 
C ^ CORRESPONDING TO J=0) . 

cr AH, H , PEL, DT , NAT 

C/; ' IT: IT= 1 KIND STRESS IS READ TO COMPUTE THE RIND STRESS CURL 

0^ IT=2 KIND STRESS CURL IS GIVEN FROM CARDS, A MAGNETIC 

C^P TAPE, OR THE SUBROUTINE CURL. 

C 

COMMON AH,DT, Y,H,REL,XM (3 9) , FJU (39) , X (3 9) , F (1 47, 3 9) ,E ( 1 47, 39) , 

1U (1 47,3 9) ,V(1 47,3 9) ,G (147 , 39) ,XU( 147, 39) , XV ( 147,3 9) ,TW{ 147, 39} ,CD, 

2BET (39) ,Y4,CCY,AR,ARY,RES (3) , DDT, 

3 MAT , L (3 9) ,LV (39) ,LA (10,37) , LB (1 0,37) , KA (10,37) , KB (10,37) , MATS , 

4MPP (1 47,3 9) ,M CQ (147, 39) ,INDIN(3) , IND(3) , NRES (3) , 

5LX1 {3} ,LX2 (3) , LY1 (3) ,LY2{3) ,IM, JM, IMP 1, JMP1 , IMH 1 
DIMENSION LI (9) ,HE(147,39) ,KN(147,39) ,C{39) 

EQUIVALENCE (U (1) ,KE( 1) ) , (V{1),WR{1)) 

READ (5,40) IM , JM, IMH, JMH 
40 FORMAT (415) 

I MP 1 =IM + 1 


IMM 1 =IM- 1 
J HP 1 =JM + 1 
CALL SHAPE1 
CALL SHAPE 1 
DO 548 J= 1 , MM 
DO 549 1= 1 , IMP1 
F (I , J) =0. 0 
G (I, J) = 0. 0 
U (I , J) =0 . 0 
- V (I, J)=0.0 
XU (1,3) =0.0 
XV (I,J) =0.0 
549 CONTINUE 
548 CONTINUE 

READ (5,18) Y, DX,D Y 
13 FORMAT (3 FI 0.0) 

P=3 . 1415926/180.0 
DX= D X*P 


MM= JM +2 

DO 100 1=1, MM 

AR= ( Y+DY*FLOAT (I) ) *P 

F 1 = C C S (AR) 

X(I)=DX*F1 
C (I) =F1 

BET (I)=F1*X (I)*14.58E-5 
100 CONTINUE 
Y=DY*P 

CD=14.58 / (DX* Y) 

DO 9997 J=1,JMP1 
XM(J) = (X (J) + X (J+1) )/2.0 
FJU (J)=CD* (X (J)-X (J + 1) ) 

9997 CONTINUE 

READ{5,16) AH, H,REL,DT, MAT 

16 FORMAT (4E15.7 ,15) 

WRIT E (6, 1 7) A H, H, REL 

17 FORK AT (1H , 1 5 REDD Y VISE 0SITY=E1 5. 7 , 1 0 X, 6 HDE PTH= E 1 5 . 7 , 1 0 X, 
115HRELAX ACCU R ACY = E 15.7) 

READ (5,10) IT 

0 FORMAT (15) 

GO TO (1,2) , IT 

1 CONTINUE 
CALL WIND 
CALL KIND 
p=0. 0 1/^- 0E5 
DO 540 J=JKH,JMP1 
X M J = X M (J) 


1 1 I *J ' I 

XM1=XM (JM1 ) 

, „ DO 5m I=IMH, III 
' IM1 =I~1 

TW{I,J)= (XMJ* (-WE (I,J)-WE (IM1 r J)) + XM 1 * (WE (I , J HI) +KB (IH1 , JH1) ) 
/; 1 + T*(WN(I # J)+WN(I r JHl)-WN(im # J)-HH(IHl # Jin)))*P 

^4 1 CONTINUE 
W TW(1,J)=TW (IKM1,J) 

Tff (2,J) “TW (IM ,J) 

TW (I MF1 , J) =T W (IMH,J) 

540 CONTINUE 
GO TO 3 

2 CONTINUE 

C CALL CURL 

CALL CURL 

R FA D { 5, 6 0 04) (TW (I,J) ,I = 3 f IM) 

6004 FOR M AT ( 5 E 1 5. 7) 

TH (IMP1 , J)=Ttf (3, J) 

TW {2 ,3) = TW (IN , J) 

TW( 1, J) =TW (IHH1, J) 

3 CONTINUE 

C INITIAL GUESS 

DO 4000 J=4,JM 

CJ= 1 * 0/ (C (J) * Y* 1 4 • 58E- 5} /4 - O 
L J= L (J) 

DO 36 C 11=1,10 
11= K A (LL, J) 

IG= KB (LL , J) 

MP= II +IG 

362 DO 361 I = II ,IG 
MM=MP-I 
K=MH+ 1 

P=F (K,J) - <TW(MM,J) + TW (K,J) ) *CJ 
F (M M , .3} = P 
361 CONTINUE 

IF (II.NE.3) GO TO 360 
F(IMP1, J)=F (3, J) 

360 CONTINUE 

F (1 , J}=F (IMMl ,U) 

F(2,J) = F (IK, J) 

4000 CONTINUE 
C . INITIAL GUESS 
CCC = 2 • 0*Y 
DO 4049 J=3,Ufl 
JP1=J+1 
LI=L V (J) 

DO 4058 LL=1,LJ 
11= LA (LL , 3) 

IG=XB (LL, J) 

DO 4053 I=II,IG 
IP1=I+1 



P= <F (I, J) +F (IP1, J) -F (I, JP1) -F (I PI ,JP1) ) /CCC 

0 (I, J)=P 
XU (I , J) =P 

P=(F ( IP 1 , J ) +F (IP1 , JP1 ) -F(I, J) -F (I,JP1) ) /XMJ 
V(I,J)=P 
XV ( I , 3) = P 
3 CONTINUE 
8 CONTINUE 

XU (I MP1 ,J)=XU (3 , J) 

XV ( IRP1 , J) =X V (3, J) 

XU(2,J)=XU(IM,J) 

XV (2 , J) = X V (IM , J) 

XD (1 , J) =XU (IM Ml ,J) 

XV (1 , J) =X V (IM M 1 , J ) 

U(IMP1,J) =U { 3 , J) 




V (1 Mi' 1 r J) = V (3 , J) 

D (2, J)=U (IM,J) 

v (2, J)=v (in, j) 

U{1 # J) = U (JMM1 ,J) 

V {1 , J)=V (IMM1 , J) 
40U9 CONTINUE 

AR= AH/6 . 37E8/6.37E8 
ARY = ?. R/Y 
CCY=2.0*Y 

V 4— Y *0-25 


RETURN 

END 

SUBROUTINE SHAPE1 

C SUBROUTINE TO BE USED FOR THE BAROTROPIC PART OF THE HOPED OCEAN 

C CIRCULATION WHEN THE SHAPE OF THE OCEAN IS CHANGED, 

C DATA:' LV, LA AND LB FOR MQQ - 

C MM, HI ( K) , M2 (K) AND M3 (K) FOR MPP. 

C IV, LA AND LB ARE THE SAME AS THOSE USED FOR THE SUBROUTINE 

C SHAPE 2 (BAROCLINIC PART), BUT THE DATA FCR MPP ARE DIFFERENT 

C FROM THOSE USED FOR THE SHAPE2 . 

C LX1, LX2 , 1Y1, LY2 WHICH SPECIFY THE AREAS INCLUDING THE SEPARATE 

C ISLANDS. 

COMMON AH r DT, Y,H, REL,XM (3 9) , FJU (39), X (39) ,F(147,39) , E (147,39) , 

1U (147,3 9) ,V(1 47,3 9) ,G (147,39) ,XU( 147, 39) , XV ( 1 47 , 3 9) , TW ( 1 47 , 39) ,CD, 

2 BET (3 9) , Y 4 ,CCY , AR , ARY, RES (3) , DDT, 

3 M AT , I (39) ,LV (39) ,XA (10,37) ,LB (1 0,37) , KA (10,37) ,KB (10,37) , HATS , 

4 MPP (147, 39) , MQQ (1 47,39) , INDIN (3) , IND (3) , NRES (3) , 

SIX 1 (3) ,I.X2 (3) , LY1 (3) , LY2 (3) ,IM, JM,IMP1,JMP1,IKM1 
I MM 1=IM- 1 
I 8P1 -IM+- 1 
J M M 1 =JM~ 1 


31 

30 


33 

34 


36 

35 


32 



JHP 1= JM+ 1 
J M H = 3 
I MH = 3 

DO 3 0 3= 1 , JM 
DO 31 1= 1 , IMP 1 
MQQ <I,J) = 2 
CONTINUE 
CONTINUE 
DO 3 2 J=3,JMH 1 
READ (5,33) M 
FORMAT (12) 

READ (5, 34) (LA (I,J) ,LB (I,J) ,1 = 1, M) 

FORMAT (7 ( 214, 2X) ) 

DO 3 5 1=1, M 
II=L A (I, J) 

IG=LB(I,J) 

DO 36 N= II , IG 
MQQ <N,J) =1 
CONTINUE 
CONTINUE 
LV (J) =M 

MQQ ( 1,3) =MQQ ( IMM 1 ,3) 

MQQ (2,3) = KQQ(IM,J) 

MQQ (IMP1 ,3)-MQQ (3 ,3) 

CONTINUE 
DO 3006 3=3,38 
READ (5, 3007) MM 
FORMAT { I 2) 

READ (5,3 008) (M 1 ( K) ,82 { K) , M3 (K) , K= 1 , M M) 
FORMAT (8 (214, 12) ) 

DO 3009 K= 1, MM 
K5=H1 (K) 

K6=M2 (K) 

K7=M 3 (K) 

DO 3010 I=K5, K6 


£f~ 


3008 



HPP1(T,J)=K7 
^3910 CONTINUE 
^3009 CONTINUE 
3Q06 CONTINUE 

DO 3011 J = 3, 3 8 
BEAD (5, 30C7) MM 
LI (J) =MM 

HEAD (5,3012) ( KA 1 (LL,3),KB1 (LL, J) , LL=1 , KM) 

3012 FORMAT (1615) 

3011 CONTINUE 

DO 6 0 03 J=1,JMP1 

MPP 1 (1, J)=MPP1(IMM1,J) 

MPPl (2, J) =MPP 1 (IH,J) 

HPP1 (IMP1 , 3) = MPPl (3, J) 

6003 CONTINUE 

DO 6004 J = 3, J M 
JP1=J+1 
U M 1 = J “ 1 

DO 6005 1=3, IM 

IF{MPP1 (I, J) . GT*1) 00 TO 6005 
IF (MPPl ( I— 1 r J) .EQ.2) GO TO 6005 
IF ( M PP1 (I , JP1 ) - EQ- 2) 00 TO 6005 
IF(MPP1 (I, JM1 ) *EQ.2) GO TO 6005 
IF(MPP1 (1+1, J) .EQ„2) GO TO 6005 
MPP 1 (I, J) =0 
6005 CONTINUE 

MPPl (IMP 1 , 0 ) = MPP 1 (3,J) 

MPP 1 (2, J) =MPP 1 (IM, J) 

MPP1 (1 , J) =MPP 1 (IM Ml , J) 

6004 CONTINUE 

WRITE (6, 6 016) (J,J=3,38) 

6016 FOR K AT ( 1 H ,5X,3 8I 3/) 

DO 6012 I=1,IMP1 

WRIT E (6, 60 1 4) I, (MPPl (I,J) ,J=3,38) 

WRITE(6,6015) (MQQ(I,J) ,J = 3,38) 

60X2 CONTINUE 

6014 FORM AT ( 1 H ,13, 2X, 3813) 

6015 FOR M AT ( 1 H ,6X,38I3) 

READ (5,2) (LX 1 (K) ,LX2 (K) , LY1 (K) ,LY2 (K) ,K=1,3) 

2 FORMAT (415) 

RETURN 

END 

SUBROUTINE WIND 

COMMON AH , DT, Y, H, REL, X M (3 9) , FJU (39) ,X (39) , F (1 47 , 3 9) ,E (147,3 9) , 

1U (147, 39) ,V(1 47,3 9) ,G (147,39) , XU (14 7, 39) , XV (147, 39) , TW ( 147,39) r CD, 
2BET (39) , Y 4,CC Y,AR ,ARY r RES (3) , DDT, 

3 MAT , L (3 9 ) ,LY (39) , L A (1 0,37) , LB (1 0,37) , KA (1 0,37) , KB (1 0 , 37 ) , MATS , 

4MPP ( 147 ,3 9) ,MQQ (1 47,39) , INDIN (3) , IND < 3) , NRES ( 3) , 

5LX1 (3) ,LX2 (3) ,LY1 (3) ,LY2(3) ,1 M, OM ,1 MP 1, J MP 1 , I MM 1 
DIMENSION WE{ 147, 7} , W N (147,7) 

RETURN 

END 

SUBROUTINE CURL 

RETURN 

END 



«S?7 
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MPP 

HPP=1 IF THE FOUR NEAREST SURROUNDING POINTS ARE WATER-POINT, 
~2 ON THE RAND (EXCLUDING THE COAST LINE) . 

FOP THE COAST- POT NTS , REFER TO FIG. 1. 
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FIG . 1 

MQQ 

HQQ = 1 IF THE FOUR NEAREST SURROUNDING POINTS ARE WATER-POINT, 
=2 ON THE LAND. 

IN OTHER CASES, REFER TO FIG. 2. 
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COMMON T (147,39,9) ,TT (137,30,5) , SP ( 14 7, 39 ) , AT (3 ) , HIT (147, 39, 9) , 

1 VV(147,39, 5) ,U (14 7,39,5), V( 147,39,9) r W( 147, 39,9) , UB (147 f 39) , 

2V8 (147,3.9) , H S ( 14 7,39) ,WN(147,39) ,Q(3?) ,XJA(3R) ,XHT(3S) , XHS ( 33) , 

3 XMY (38) , % ( r ) ,X (39) , XM (39) , FU(39) ,?jn (39) ,77 (6) ,J], A7, A H , DT , AT., 

4G AM , Y,YY , R , PR , HH, H2,AZHH, R2, A BY , AHRF, AH2, B3 , BE? , A A A , EBP , DDT, A 1 , UV, 
5UUIT , V VV , R R H , A 0 , G A M 2 , A L 2 ,B UN , E Y, AM , A MR R , GA EGA M, AT AL, BAN , KMP1, 

63 PP ( 1 47, 3 9) ,KQQ (1 47,3 9) , I 3 , 33 , K S* , I M Ml ,3 EMI , KM FI , IMP 1 , J F P 1 , M AT, 

7 NK, N NN N , I HH , JMH , MATS ,L (39) ,KA ( 10,39) , KB (10, 39) , LV (3D) , HIND, 

8LA ( 10,39) ,T.B (10,3 9) ,NN , TE (14 7,46) ,UN( 147, 39) , TO ( 147,46) ,PM( 147, 39) 
* , CL ( 147, 3 9) 

DI ME NS I C N UH E AD ( 3 9) , TH EA D (3 9) 

DIMENSION QF (147, 39) , QH (1 47,39) ,QSH (1 47, 39) 

CCMMCN /TEITP/ Q F , Q H , 0 S H , 2 0 E , 7 QH , 7.S H , 7 TR N 
DIMENSION FOE (39) ,7QH(3 9) 

LOGICAL STOP 


PEAL KK 

DIMENSION XRN (39) ,I$H (39) 

DIMENSION IT (39,9) ,79 (39,9) , TJ Z ( .3 9 , 5) , VZ ( 39 , 9) 

EQUIVALENCE ("1(1) ,U7,( 1) ) , (ZW { 1) ,VZ ( 1) ) , (I?( 1) ,9(1)) , (IV (1) , W (19 6) ) 
DIMENSION ICTL(7) ,P.CTL (2) , FIELD (18) 

DATA IC T L/ 1,6, 9,1, 18,7,2/ 

DATA RCTL /0. 001,60./ 



fill i' UK ! y| A i. ( 1 ) ’ , i 'J A , ’ At’hliii'i r.» L -> i C a l i.JJ ) 

UHEAD (1) ~ - 17 . 

THEAD(I) = -74. 
v ^ DO 88 1=2 ,39 

^ * UHEAD (I) = U If F* A D ( 1 ) + ft .*(1-1) 

*-\ e THE AD (X) = TUFAP{1) + 4.*(T-1) 

m CONTINUE 


NNNN=0 

DEAD (5, 9) KT 
READ (5,9) TN 
9 FORM AT (15) 

IF(IN.EQ.O) 00 TO 11 
17 READ (KT, FNU=92) 

* TT , UU , 7 V , T , a , V, .AT , S P, W , US , V R, a, X H r Q , FU , FI TT , P? , A Z , AH , 7 F , 

1 W W , B A N f A L , G A H , Y , Y Y , E , fi R , H H , H 2 , A Z R H , F [1 K , R 2 , A HY , A H P P , A K 2 , P B , F 3 2 , 

2 AAA, 3SB, DDT, A 1 , 117 , 11 UU , V VV , R RH , A Q, GAH2 , AL2 , SY , AM RP , GAHGA M , ALAI, 

311 ,Z , Z 7, A M , X J A , X IIT , XHS , X MY , CD, 

UMPP,MQQ,KA,K3,L,LA,LB,LV,MATS # NFND r NK,tfNNN,8AT,IH,JM,KM,IMH f r M 3 1 , 
5 T M P 1 , K M P 1 , K M M 1 , J M M 1 , J HIM , N N , J H H , M A 0 , T E , T D , R W 
6,QE,Q.H,05K,CL,UN 
GO TO 17 

92 BACK SPACE KT 

IF (MAG. LT .34) GO TO 94 
REWIND KT 
KT = KT + 1 
HAG = 0 

94 CONTINUE 

7 3 READ (10, EM D-7 4) OS,VB 
GO TO 73 
74 REWIND 10 


DO 2222 1=1, 1 45 

WRITE (5,7000) I, (TIN (I ,J) , J=3,20) 
2222 CONTINUE 


DO 4 4 44 1 = 1, 1 45 

WRITE (6,7000) I, (UN (I ,J) , J=21,38) 
4444 CONTINUE 

WRITE (6, R 9 0) (THE AD (J) ,0=3,2 0) 

DO 7040 T=1, 1 47 

WRITE (6 ,7000) I, (TE ( I , J) , J= 3, 2 0) 
7040 CONTINUE 

WRITE (6, 890) (THE AD (J) ,J=21,3R) 

DO 7034 1=1,147 

WRITE(6, 7000) T, (TE(I,J) ,J=2 1, 3B) 
7034 CONTINUE 

CL (1 ,3) =6 .00 
CL (2, 3) =6. 00 
CL ( 1 4 5, 3) =fi..0 0 
CL (1 46,3) = 6 . 

CL ( 147, 3) = 6. 

DO 3333 T=l f 1 45 

WRITE (6,7 000) I , (CL (T , J) , J= 3, 20 ) 
3333 CONTINUE 

DO 5555 1=1,145 

WRITE (6,7 0 00) I, (CL(I,J) ,J=21,38) 
5555 CONTINUE 

FA = 4 00000 . 

DO 311 1=1,147 
WE (I, 3) =1 . 16* FA 
WE (I, 4) =1.72* FA 
WE (1,5) =1 „ 8 3 * F A 
WE (1,6) =1 .90*1 A 
W E ( T , 7} =1 . 8 3 * F A 
WE (I ,8) =1 .50* FA 
WE (1,9) =1 . 10* FA 
WE (I, 10) =0.7 2* FA 



w { J , I I ) — U B ■ > U J. i ! 

WE (I, 12) =~0.07* FA 
WF (I, 1.3) =-0. 7 0* FA 
v * WE (I, 14) =-0.5 2* FA, 

*' " WE (1,15) =-0.6 0*FA 
WE (I, 16) =-0.00* FA 
W E ( I , 17) =-0. 5 2*?A 
WE (1,18) =-0.3 8* FA 
WB (I, 19) =-0.30* FA 
WE( I, 20) =-0.2 4* FA 
WR(I,21) =-0. 2 6* r A 
WE (1 , 22) =-0. 3R*Ffl 
WE (I, 23) =-0. 9 0* FA 
WF { I , 24) =-0.5 3* FA 
WE (1,25) =-0.4 0 * P A 
WE (1,26) =-0.1 4* FA 
W E { I , 27) =0. 10*FA 
WE (1,28) =0.42 *FA 
WE (1,29) =0.61* F A 
W E ( T , 30 ) =0.7 4* FA 
WE( J , 31) =0.78*'^A 
WE (1,32) -0.63*? A 
WE (I, 33) =0.49* FA 
WE (I, 3 4) = 0. 14*PA 
WE (1 , 35) =+0.01 *FA 
WE ( I , 36 ) =- 0.04* FA 
WE( I, 33) =- 0 - 1 3* FA 
311 CONTINUE 
A H= 2 „ 5.F 7 
A H Y = A H * Y 
AHRR = AH/RP 
A H 2 = A H/ ( Y Y* R R ) 

AL= 0 . 25 E- 3 
ALAL=AL*2.0 
AL2=AI/2. 0 
A E= 1 . 0 
AZ H U = AT / H H 

" ' GANGAN=2 5.0*H/{AZ*0.864 E5 ) 

WRITE (6 , 7777) 

WRIT E ( 6, 9006) 

9005 FORM AT { 1H0 ,21 H3AHCTR0PIC COMPONENTS/) 
WRITE (6,9110) 

.9110 FORM AT ( IX, 1 03’ ) 

WRITE (6 , 9004) (J,J=3,21) 

9004 FORMAT (1H0,4X , 19 (2X,T2, 2X) /) 

T)0 9 0 00 1=1,1 47 

WRTTE (6,9001) I , ( U B (1,0) ,0 = 3,21) , I 

9000 CONTINUE 
WRIT E ( 6 , 7777) 

WRITE (6, 9004) (I , J =2 0 , 3 3 ) 

DO 9002 1 = 1 f 1 47 

WRITE (6, 9001) I, (HB (I , J) ,0=20,38) ,1 
9002 CONTINUE 

WRITE (6 ,7777) 

WRITE (6, 9 200) 

9200 FORMAT ( IX, * VB 1 ) 

WRITE (6 , 9004) (0, 3=3,2 1) 

DO 9006 T= 1 r 1 47 

WRITE (6 , 900 1) I, (VB (I r J) ,0=3,21) ,T 

9006 CONTINUE 
WRITE (6,7777) 

WRITE (6, 9004) (.1,0=2 0, 38) 

DO 9007 1= 1, 1 47 

WRITE (6,9001) I, (VB (I , J) , J=?0,39) ,T 

9007 CONTINUE 

9001 FOR M AT ( 1 H , 13 , 1 Y , 1 9 F6 . 2 , 1 X , I 3) 



Vh< L'i h (b , 'MDU) UM JJ ,^-3,-iH) 

9100 .FOR M AT ( 1 H 0 , 1 9 F6 . 7/1 H ,18F6.2/) 

IF(NK.NE.l) GO TO 4 92 
v V GO TO 5 
^,11 CONTINUE 
*5 N N = 0 

CP= 0 « 0 
NK= 3 
M A G = 0 
CALL FREP 
CALL FRF PA 
CALL GEO ST 

5 READ (5, 3) EAT S, REND, NPP.XNT, NT A PE, DT 

3 FORMAT (4 IS, SI 5.7) 

WRITE (6, 6) N PRINT , NT A. PE 

6 FORM AT ( 1 NPPINT = *,15,* NTAPB = * ,15.) 

WRITE (6,9) KT , HAG, CD 

8 FORMAT (* KT = r ,T2, f HA G = *,1?,* CD = * , F7 . 1 ) 

DTD=DT/0. 8E4.ES 
K = N N + 1 
N EN D— NEN R + K 

WRITE (6, 4) K, MEND, HATS, DT , A 2 , AH, A H 

4 FO P H A T ( 1 1 , * T THE STEP FROM ',15,37,^0 * , IS, 5T , 1 M A TS= 1 , 1 3, 
X * DT= f , UE IS. 7/) 

IF(t>T.LT. 10.) STOP 100 

DDT=DT*2 . 0 

DO 1 1 1 J = J H rf , J H M 1 

X J= F JU (.1) *DDT 

XJA ( J)-XJ*AAA 

XHT (J) = 1 . 0/ (1 .O+XJA (J) **2) 

XHS (J)=XJ*BBB 
XHY ( J ) - V * X H (J) 

111 CONTINUE 

IF (NN.GT.O) GO TO 251 

CALL CONTI 

MAT-2 

WRITE (2) TT 
REWIND 2 
WRITE (3) U,V 
REWIND 3 
DDT=DT 
GO TO 253 
"251 NN=NN+1 

CD=CD+DTD 

HHN“ (NN/HATS) *EATS 
IF (NNN. NS. SN) GO TO 2 SO 
M A?= 1 
WRITE (6,2) 

2 FORMAT (1 H ,5H BEGIN) 

2 SO CONTINUE 
253 CALI, THERMS 
CALI, ,£Q A 
CALL CONTI 

WRITE (6,7013) THJTJ , V V V, TIV , UN , CD 

7013 F OEM AT ( 1 H , 7 IT ENER G Y=3 {E 1 5 . 7 , 3 X) , 5HSTE P-15 ,2X, 411 DAY-F 10 . 2) 

IF ( KA7.NE. 0) GO TO 70 19 
CALL WTO (STOP) 

IF (STOP) N FN D = 1 0 

.TF ( NN.GE. N.END) GO TO 70 15 

IF ( ( (NN/NTA^E )*HTAPE J.NE.NN) GO TO 7025 
IF(NN.EO.O) GO TO 7 019 
7015 HAG = HAG+1 
WTM = 0.0 
GQSH =0.0 
GQE — 0 .0 



b K N ~ U . 

DO 20 C 1 J = 3 , 3 8 
** DXYP = X(.1) * Y 

CL AT = A 5. S (DXYP) 

^ KE=0. 

ZQH (J) = 0. 

ZQE (J) = 0,0 
ZZZ=0. 

7 RN (J) - 0. 

DO 2 0 02 1 = 3,1 46 

IF ( N PP (1,0) - GT . 1 . A ND . HPP (I, J) . I.E . 5 ) GO TO 0011 
UM=UN (I, J) *51 .48 
TE I J=?E (T t J) 

TIJ= T ( I , J , 1 ) 

TW = TD (I, J) + 2? 3. 15 
12 3 TS=T (I , J , 1 ) +2 73 . 1 5 
PCB=1013. 25*0.1 
ESS = 1 0 . 0 * * ( 8 . 4051-2353. 0/TS) 

ESD = 10. 0** (8. 4 051 -2 35 3 . 0/TW) 

QS=0 . 622* ESS/ (PCB-ESS) 

QA-0 . 622*ESP/ {^CB— ESD) 

IF ( ( (TIJ-TFT. J) +0, 609* TS* (QS-QA) ) . LT. 0. 0) GO TO 33 4 
332 CU= 0.033 

CO= 1 . 23*0. 001 
GO TO 133 
3 34 01.1 = 0 „ 01 8 

CO=0 - 41*0. 001 

1 33 Qfi (I ,0) =- 1.2*0.00 1 * . 2 4 * CO *U»* (TIJ-TEIJ) *0.864 F5 
OE (I, J) =-1 .2*0. 00 1 *CO *595 . *U 11* ( OS -QA) *0.864F5 
RN (I , J) = RN (I, J) *0 . 864 F5 
OSH (I,J) = RN (T ,,1) + OH (I, J) + QF (I, J) 

GQSH = GOSH+QSH (I , J) *CLAT 
GQE = GQE + OF (I, 3) *CLAT 
GOH = GOFi -+ Q H (T , .3) *C I. A T 
GEN = GRN + RN(I,J) * CL AT 
STH = KTfl+CLAT 
KK = KK+ 1 . 

ZZZ = QSH (T , J) + 777 

ZQH(J) = Qfi(I,J) + ZOH(J) 

ZQE(J) = QE(I,J) + ZQF(J) 

ZRN(J) = RN(T,J) + 7RN(<1) 

GO TO 2002 
9 011 QSH (1,3) = 0. 

QE ( I , J) =0. 

Qfi (1,0) '0. 

2002 CON TIN (7 F 
KK = 1./KK 
2SH ( 3) = ZZZ* KK 

7QE(J) = ZCF(O) * KK" 

ZOH (J) = 7. OK ( J ) * K K 
ZRN (J) = ZRN(J) * KK 
2001 CONTINUE 

NTH = 1./WTH 
GQSH = GQSH* WTZ 
GQE = GOE*VfTM 
GQH = GQH* WT 3 
CRN = GR N* WTM 



* ,QE, QH, QSH, CL, ON 
FNDTT'TT.F. KT 


1 

h l. 

^ *- 

7*02 5 

W! 4 
7020 


20 3 

301 

300 

7010 

7012 

7011 


7022 
■7 009 

7053 

7021 

390 



BACK SPALL* aV 
WRITE (6,1) MAG 

FORMAT ( 1 FT , 3 5 HNUMRE R OF RE CO? OS IN MAGNETIC TAPE=I2) 
IF (NN.GE. NEND) GO TO 7014 
IF (MAG. GF .34) GO TO 93 
I F ( TJ V . GT . 1 . 0 E 2 6) GO TO 7014 

IF ( ( (NN/N PRINT) PRINT) .HE* NN) GOTO 7019 
WRITE (6,7 020) NN, CD 

FORK AT (1H ,1 1 IiTTMS (STEP) =18, 5X, lOHTIME(DAY) =E15.7) 

DO 300 K = 1 , 5 
KK= 0 . 

DO 301 ,7=3,33 
UZSUM=0. 

V7.SU M=0. 

DO 3 03 1 = 3,146 

I F ( M QQ ( I , J) , I> E « 5 . A HD . II Q 0 ( I , J ) - N E . 1 ) GO TO 30 3 
KK-KK+1. 

U7S U M=U ZS U M + U (I , J , K) 

VZ5 fl M=V 7.S fJU+ V (T,J,K) 

CONTINUE 

IF (KK .IT. 1 .) GO TO 3 00 

II 7 (J, X) =U7Sim/KK 

V7 (J , K) = VZSDM /KK 

CONTINUE 

CONTINUE 

DO 7 0C9 K = 1 , K M 

WRITE (6, 7777) 

WRITE (6,7010) K , CD 

FORMAT (1 X, MJU AT K = * ,11,* DAY *,F7.1) 

WRITE (6,890) (UILEAD (J) , ,3=3, 20) 

DO 7011 I = I M H , ID 

WRITE (6,7012) I , (OU (I , J,K) ,J=3, 20) 

EOF M AT ( 1 1I ,13 ,1X, 18F7.2) 

CONTINUE 

WRITE (6,7778) ( T J 7 ( J , K) , ,7=3 ,20) 

WRITE (6,7777) 

WBITE(6, 7010) K , CD 

WRITE (6,8 90) (UHEAD(J) , ,7=21,38) 

DO 7 022 T = I E H , T E 

WRITE (6,7 012) T, (UTJ (I, J ,K) , J=21 ,3 8) 

CONTINUE 

WRITE (6,7778) (U7( J,K) , J=21 ,3 0) 

CONTINUE 

DO 7 0 70 K = 1 , K M 

WRITE (6,7777) 

WRITE (6 ,7 053) K , C D 

WRITE (6, 9 90) (TJHEAD (J) ,,7=3,2 0) 

FORM AT ( IX, 1 TV AT K = * , 1 1 , ' DA Y * , F7 . 1 ) 

DO 7021 T = T MH , T M 

WRITE (6,7012) T, (VV(T,J,K) ,,1=3,70) 

CONTINUE 

WRITE (6,7 778) (VZ ( ,7 , K) , J=3 ,20) 

FOR MAT (/5X, 13F7.2/) 

WRITE (6, 7777) 

WRITE (6,7 053) K,CD 

WRITE (6 ,890) (UHEAD (J) , ,7=21 ,33) 

DO 7023 T=imi,IM 

WRITE (6,7012) I, (VV(I,J,K) ,0=21,38) 

CONTINUE 

WRITE (6, 7778) (VZ(J,K) ,J = 21,38) 

CONTINUE 
DO 200 K=1 ,KM 
DO 201 ,7 = 3,33 
KX = 0 „ 

WZSUM=0 . 

TZSUM=0 . 


7070 



DO 2U3 1 = 3, {4b , 

IF(KPP(I,J) * R T . 1 * AND. HPP(J, J) . LE. 5) GO TO 203 
w KK=KK+1. 

V WZSUM=WZSUH+W (T,J,K) 

^ r TZ ST1 M=T ZS OK+TT (T , J , K) 

3 CONTINUE 

, IF ( KK.J.T. 1 .} CO TO 200 
ZH ( J ,K) = WZSn?1/KK 
ZT (0 , K) =TZSUM/KK 
201 CONTINUE 
200 CONTINUE 

DO 7 OOF K = 1 , K M 
WRITE (6 ,7 777) 

WRITE (6,7006) K 

7006 FOFMAT(1K , 1 7 HTEMPERATURE AT K=I1) 

WHITE (6, BOO) (THE AD (J) , J=3,20) 

DO 7 050 I = I HR , I M 

WRIT E (6,7000) I, <TT <1 , J ,K) ,J=3,20) 

7050 CONTINUE 

WHITE (6,7778) (ZT(O r K) ,J=3,20) 

WRITE (6 ,7777) 

WRTTE(6, 7006) K 

WRITE (6 , 890) (TREAD (.1) , J=2 1 ,38) 

DO 7 051 I = I H H , I M 

WRIT E (6,7000) I, (TT (I, J ,K) ,J=21 ,3 8) 

7051 CONTINUE 

WRITE (6, 7773) (ZT(J,K) ,J=21,38) 

7005 CONTINUE 

DO 7 0 02 K~2 , K M 
FK = K-.5 

WRITE (6, 7777) 

WRITE (6, 7001) EK , CD 

7001 FORMAT (1 7, * VERTICAL VELOCITY AT K = ’,F4.1,« DAY ’ , E7 . 1 ) 
WRITE (6,890) (TREAD (J) , J=3,20) 

DO 7003 

WRITE (6 ,7000) I, (W (I, J, K) , J=3, 20) 

700 0 FOR EAT (If! , 1 3 , 1 Y , 1 8 F7 . 2 ) 

7003 CONTINUE 

WRITE (6 ,7778) (ZW(J,K) ,J=3,20) 

WRITE (6,7 777) 

WRITE (6, 7001) FK , CD 

WRITE (6, 890) (TREAD (J) ,1=2 1,38) 

DO 7 024 T-3Sfi,IR 

WRITE (6, 7000) I, (W(T.,J,K) ,1=21,38) 

7024 CONTINUE 

WRITE (6,7778) (7W (J, K) , J=21 ,38) 

7778 FORMAT (/3 0X, 1 ZONAL MEANS 1 //5X , 1 8F7 . 2) 

7002 CONTINUE 

IF (UV.OT.1.0E26) STOP 7 92 
WRITE (6,7777) 

WRITE(6, 9016) CD 

9016 FORMAT ( 1 X, 1 SURFACE HEAT FLUX (CAL PER DAY) DAY = » ,F7.1/) 
WRITE (6,9020) GQSH 

9 020 FORMAT (* G TOR AL MEAN = • , FI 0,3) 

CALL PRINT (QS H , 7-S H , TH E A D , ICTL, RCT.L) 

WRITE (6, 7777) 

WRITE (6 , 9017) 

FORMAT {’ LATENT HEAT FLUX’) 

™ WRJ TE ( 6, 9 020) GQE 

CALL IRI NT (QS , ZQE ,T HEAD ,1 CT L, 9CT L) 

WRITE (6,7777) 

WRITE (6, 90 18) 

9018 FORM AT ( ’ SENSIBLE HEAT FLUX’) 

WRITE (6, 9 020) GOH 

CALL PRINT (OH, ZQH ,T HEAD , ICT L, POT L) 

WPT? F (6 , 7777) 


<>023 

V 4- 
♦ t 

t h- 

v7i)1 9 
7017 

20 

21 

22 

250 1 
89 

X 

X 

3 005' 
3004 

3001 

3006 

3007 

3002 
3 008 
3 009 

m 


W H.L lh {h e 9U2 J) 

FORMAT** SURFACE RADIATION FLO X (CAL/CM**2/P A'{) *} 
WRITE (6,9 020) GEN 

CALL PRINT (RN , TEN , TREAD ,ICTI , RCT L) 

IF (NN.GE. NEED) GO TO 2501 
CONTINUE 

IF ( MAT. EO. 0) GO TO 25 1 
MAT = MA? + 1 

IF ( tf AT. LT . 8) GO TO 7017 
DDT=DT*2. 0 
GO TO 251 


DDT -OT 

IF( MAT. EQ.3) GO TO 20 
IF (MAT. EO. 4) NK=3 
IF (MAT. EO .5) GO TO 21 
IF ( MAT. EQ . 6) »K=3 

GO TO 250 
NK= 1 
MAT' =7 


NK= 3 

GO TO 22 
NK = 2 


NN=NN+1 
CE-CD+DTD 
GO TO 250 


CONTINUE 

DO 3 0 C5 J -3,37 

LJ = L V (J) 

DO 3005 r 1 1 = 1 , T.J 
II = LA (Ml, .7) 

FORMAT (* 0* ,nx , 19F6. 1/) 

TG = LR(M1,J) 

DO 3005 I=II , JG 

UtT (1,3, 1 ) - U U (T, J, 1) *? s (1 ) +DU (I ,3,2) *7 (2) + UU (1,0,3) *7 (3) 
+ U!J ( J, J,4) *2 (4) + U5J (1,3,5) * 7> (5) 
w (I, J, 1 ) = VV(I , J, 1) *Z (1) +VV (T,J,2) *Z (2) + 7V (I, J, 3) *T (3) 
♦ VV (I ,.7,4) *7. (4) + VV (I, J,5)*Z (5) 

CONTINUE 

FORMAT* 1 * ,13 , IX, 19F6. 3 , 1 X,I3) 

WRITE (6,7777) 

WRITE (6, 3001) 

FORMAT* * UUBAI? » /) 

WRITE (6,89) (UKFAD (J) , J=3,21) 

DO 3006 1=3,146 

WRITE (6, 3 004) I, ( U U ( I , J , 1 ) ,J=3,21) , I 
WRITE (6,0 9) (U FIE AD (J) , J=-2 0, 3 8) 

DO 3007 1=3,146 

WRITE (6,3 004.) T, (UU (I, J ,1 ) , J=20 ,3 0) ,1 
WRITE (6,7777) 

WRITE (6, 3002) 

FORMAT {* WEAR*/) 

WRI T E (6 , 3 9) (AHEAD (J) ,.7=3,21) 

DO 3 0 OB 1 = 3, 1 46 

WRITE (6, 3 004) I , ( VV (I , J , 1 ) , J= 3, 2 1 ) , I 
WRITE (6,8 9) (IJHFAD(J) , J = 2 0,3B) 

DO 3009 1=3,146 

WRITE (6,3 004) I, (VV(T ,J,1) ,J=20,3 8) ,1 
STOP 


RFWJUD KT 
KT=KT + 1 
MAG = 0 
GO TO 7015 


END 

SUBROUTINE PRINT ( FI, 7, S P, T UFA 0, 1 CT L, RCTL) 

DIMENSION THE AD (3 9) 

DIMENSION FIELD (18) ,7S? (3 9) , p I ( 147,3 9) ,ICTL(7) , RCTT ( 2) 



itn i URn A .1. ( ' I 1 , 3U A , * r .VO. Hi m2, x l jU 'j 
7012 FOE M AT { 1 1 , 1 3 , 1 X , 1 8F7 . 2) 

3 90 FORMAT {/5X,1 877*2/) 

7778 FORM AT (/30X, ' ZONAL ME A NS » //5X r 1 8F7 . 2) 
,H ’ WRITE (6 , 890) (THE AD (J) , 3=3,20) 

CALL SHADE (ICTI, PCTL, 0 . ) 

DO 9012 3, 1 46 

DO 701 J - 3 , 2 0 
3 1 = J- 2 

tv I FIELD (J1) - 71(1, J) 

CALL SHADE (0,0. r FlELD (1) ) 

WRITE (6, 7 012} I , FIELD 
9013 CONTINUE 

WRITE (6 r 7778} (ZSP(J) ,0 = 3,20) 

WRITE ( 6, 7777) 

WRITE (6 , 890) (THEAD(J) ,J=21,38) 

CALL SHADE (I CTL, RCTL, 0 . ) 

DO 9015 1=3,146 
DO 702 3=21 ,38 

3 1 = 3-2 0 

702 FIELD (3 1 ) = *1(1,3) 

CALL SHADE (0, 0* , FIELD (1 )) 

WRITE (6,7012) I, FIELD 
9015 CONTINUE 

WRITE (6, 7778) (ESP (J) , 3=21 ,33) 

RETURN 



M A R = 0 

40 IF (KP1 -EQ.KM) GO TO 50 
K = KP1 
KP 1=K+1 

B= (AT (K) ~AT (KP1) ) / (ABS (AT (K) ) ) 

IF (B .GE. -5, OB-6) GO TO 40 
D Z 1= Z (K) 

D7.2 = Z (KP 1) 

A VT= (AT (KP 1 ) *DZ2 + AT (K) * DZ 1 ) / (DZ2 + DZ1) 
AT (K Pi) = A VT 
AT (K) =A7T 
MAR = 1 


GO TO 40 

5 0 IF(HAR.FO.O) GO TO 51 
KR=MR*1 

I F ( M R . L T . 1000) GO TO 41 
WRITE (6,1) (AT ( T ) ,1 = 1, KM) 

|1 FORMAT ( 1 H ,23 HTOO MANY ODER IN ADJUST , 5E15*7) 

T>1 RETURN 
END 

SUBROUTINE EO A 

COMMON T (1 47, 3 9,5) , TT (1 47,39,5) ,SP( 147,39) , AT (5) ,110 (147,39, 5) , 

1VV ( 1 47, 39, 5) ,3( 147,39,5) ,V (147, 39,5) ,*(147,39,5) , UP (147,39) , 

2 VB ( 1 47, 3 9) ,9F (147 ,39) ,WN( 147,39) ,0(3°), X J A { ? 8 ) , XHT (38) , XPS ( 38) , 

3 XMY (38) r Z (5 } , X ( 39 ) , XM (3 9 ) , PU ( 39 ) , FJU ( 39) , ZZ (6) , H , A? , A II , DT , A L, 

4 GA M - V , X Y ,R,RR!,HH, H 2 , A Z Hll , F 2 , A H 7 , A H P P. , A 82 , r 'B , B P9 , AAA , ?RB , DDT , A1 , UV, 


)\ ' J I 


f f 

V t 


••> UUU , V VV , UHH, AQ, GA wl , AJ.2 ,.u>.l h, ri A , iiri , An Mi , o *\ i it-. A . !, A i.a L, ».•« I, 

6MPP (1 47, IQ) , MOO (1 47 f 3 0} , T fl , .7!‘: , K 5, IP Ml , JM tt 1 , KM Ml , IMP 1 , JFT>1 , M AT, 

7 NX, N NN N , I M f I , J M II , MATS,!. (30) , K A ( 10,3°) , KD (1 0,3 9} ,LV{^0) , NFNP , 

81. A (10,30) , 173 (10,39) , N N , TE (147,46) ,0N ( 147,39) , TJ) (147,46) , PH (147, 39) 
* , CI, (147,3 9) 

DIMENSION SEN (147,3 9,5) 

EQUIVALENCE ( W ( 1 ) , SEN ( 1 ) ) 

C CM MON/TEM P/T FN (146,18, 5} 

DO 3000 J= JMU , JMM1 
JP1 -.7 + 1 
J M 1 = J - 1 


X M <_7 - X M ( J) 

X M J 1 = X M ( J Ml ) 

XJ=X (J) 

XJP1=X (J?1) 

XJ1 = X (3 Ml) 

AD -A 1 *XJ/XBJ*0.3 
AC=A1*XJP1/XMJ*0. 3 
A B= A MHR/ { X S J* X M J) *0.3 
R XM J=R*X« J 


L<1~LV (J) 

DO 4003 H 1 = 1 , L J 
11= J, A (Ml , J) 

IG= LB ('ll , J) 

DO 3001 1=1 I, 10 
MQ= MQQ (I , .7) 

IP1 =1+1 
IM 1=1-1 

C*****cOF!PT.mi?CN OF PRESSURE 
GO TO (1,2,3), NX 

3 CONTINUE 

SP0 = SP(I,J) -SP(IP1 , JP1) 

SP1 = SP (IP 1 , J) -3P(I,J?1) 

PXB IJ=3P 1— SPO 
PYBI J=-SP1 -SPO 
GO TO 4 

1 PXB T.7=(SP(IP1 ,JP1) -SP (T,JP 1) ) *2.0 
PYBIJ= (SP (IP1 , JP1 ) -SP (T Pi , J) ) *2.0 
GO TO 4 

2 PXB I v.7 = ( S F ( IP1 ,J) -SP(I,J)) *2.0 
PYBIJ= (SP (I,. 7 PI ) -SP (T ,J) ) *2.0 

4 CONTINUE 

C-*****cOMP!)TATTOH 0^ PRESSURE 
WEI J=WE (I, J) 

WNT J=WN (I, J) 

DO 3 002 K=1,KH 
KP 1 =K + 1 
2 K= 7, ( K) 

2Z1 =ZZ (KP1 ) 

KH1=K-1 
Z7.K-ZZ (K) 

5555 HS= V (I, J , K) 

HT= U ( T , J , K) 

H Ml J = U (I Ml , J, K) 

VM1 J = V (I'll,, 7, K) 

UP1 0 = TJ ( T P 1 , J, K) 

VP 1 J = V (IP1 ,J, K) 

UIP 1=1) (I, JP1, K) 

VIP 1=V(I , OP 1 , K) 

DIM 1 = U (J , JK1 , K) 

VIM 1 = V { J , J M 1 , K) 

5556 UIJ = rTT 
VIJ=HS 

AX=AZHH/ZK 

IF (K . EO. 1) GO TO 3003 
DII7.= (V f T , .7 , K M 1 ) -HT) /7 7.K 


\) V h - ( V ( 1 , .J f t\ y -i \)~l- o) / L 
IF(K.NE.KM) HO TO 300 U 

► - duzp 1 = 0.0 

dvzp i-o . o 

^ 1 GO TO 3 0 09 
*3003 DOZ = WEIJ 
A 0 VZ- WNI J 

DUZP 1= (HT-0(T # J # KP1) ) /ZZ1 
DVZP 1= (HS- V(T , J,KP1) ) /ZZ1 
3 009 HT=?iX* (DOZ-DUZPI-WFIJ^ZK) 
HS= A X* ( D V Z- DV ZP 1 - W HI J*Z K) 


IF(KQ.EQ. 1) GO 

TO 

300 8 

IF (MQ.EQ. 1211) 

GO 

TO 

12 11 

IF (MQ.FQ. 1 112) 

GO 

TO 

11 12 

IF (MQ.FQ. 1121) 

GO 

TO 

11 21 

IF { KQ. EQ . 2211) 

GO 

TO 

2211 

IF( MQ.FQ. 2 112) 

GO 

TO 

21 12 

IF (MQ.EQ. 1221 ) 

GO 

TO 

1221 

IF (MQ.FQ. 1 122) 

GO 

TO 

1 1 22 

IF(MQ.EQ. 2121) 

GO 

TO 

2121 

IF(MQ.EQ. 1212) 

GO 

TO 

12 12 

IF (MQ.FO. 2221 ) 

GO 

TO 

2221 

IF (MO.FQ . 221 2) 

GO- 

TO 

2212 

IF (MQ. EQ. 1222) 

GO 

TO 

1222 

IF (MQ. EQ. 2 122) 

GO 

TO 

2122 


2111 UP1J=-0TJ 
VP1 J=-VIJ 
GO TO 3000 
1211 TJIP1=“UIJ 
VIP 1=-VI J 
GO TO 3000 
1121 UM1J=-TJU 

vm j=-vij 

GO TO 3008 
1112 VIM1=-VIJ 
UIM1=-UTJ 
GO TO 3008 
2221 TJP1J=-UTJ 
VP1J=-VIJ 

uipi=-hi j 

VIP 1--VIJ 
UM1 J=-UIJ 
VH10--VU 
GO TO 3008 
2212 UPl.7-=-UTJ 

vpu=-vtj 

VIP 1 — VIJ 

mpi=-ni j 
uih i=-ni.T 

VIM 1 =-VI J 
GO TO 3008 
122 2 UIP 1 =-IJI J 
VIP 1 = -V I J 
. 0 fll *1 UI J 
V Ml J -- V I J 
0IM1=-UIJ 
VIM 1 --V.I J 
GO TO 3008 
r 22 U P 1 J =—01 J 
VP1 J = -V T J 
tJWl J--UI.7 
VH1 0=-VI.l 
UIH 1=-UIJ 
VIM 1 =-VI J 
GO TO 3008 
2211 nPU=-UTJ 



V Pi J -~V JL.J 
UIP1--UIJ 

^ VTP1=-VIJ 
GO TO .3 0 08 
‘^321 UP1 J = “UI J 

V PI J = -VT 3 
UH1J=-0TJ 
VK1.7--VIJ 
GO TO 3 00 3 

2112 UP1.1=-UIJ 

V P 1J = - V I J 

uim=-ci j 

. VIM1--VI.1. 

GO TO 3 008 
1221 U IP 1 =- U I J 
VIP1=-VXJ 

unu=~uij 

VM1 J=-VJJ 
GO TO 3003 
1212 UIP1=-TJT.1 

VIP1=-VIJ . 

UIM 1 =~UT J 
VIM 1 =— V T J 
GO TO 3008 
112 2 UMIJ^-rilJ 
VK1 J=-VI J 
UIfH=-UI.T 
VI M 1 = -VT 0 

3008 HT-ITT+AB - ( UP 1 ,1 + UM J-2 . 0*OIJ) + AC* (HIPI-HU) - At) * ( UT J-M TH1 ) 

HS=HS+AB* (VP 1 J+ VN 1 J-2 . 0*VIJ) +AC* (VTP1 -VIJ) -A D* ( V J J- VI Ml ) 

IF(K , SE. 1 ) GO TO 8000 
P A = 0 . 0 
PB=0.0 
GO TO 8001 

4000 GO TO (5,6,7) , N K 

5 PA=PA+ZZK* ALAI*{T (1,001, K ) +T (I r JPl,Kf!l ) -T (I PI, OP 1 , K) -T (T P 1 ,,1P1 , K 

1M1) ) 

P3-PB+ZZ K* ALAI* (T {TP 1 , J,K) +T (I Pi , J, K H 1 ) -T ( IP 1 , .IP 1 , KJ -T (I PI ,8 Pi , K ’1 

11)5 

GO TO 4001 

6 PA=PA+ZZK* ARAL* (T(T, J , K) +T (I,J,KF1) -T(IP1, J , K) -T (TP 1, J ,KF!1) ) 
PB=--P8+ZZK* A L A L * ( T (I, J , K) + T (I , J , K HI } - T (I , J? 1 , K) -T (T , JP 1 , KM 1 ) 1 
GO TO 4 001 

7 CONTINUE 

P 1- T ( T. , 0 P 1 , K) +T (I,JP1 , K Ml ) -T {IP1, J,K) -T (IP1, J,KK1) 

P2=T (T,J,K) +T (I, J, Kill) -T (T?1, JP1, K) -T (IP1, JF1,KH1) 

P A = P A +Z 3 K * AL* (P1 + P2) 

PB = PB +7, 3 K * AL* (P2-P1 ) 

4001 TEN (I,J,K)= < (-FA* AQ>PXBIJ) / R X BJ + HT ) * DDT 
SEN (I,J, K)= ( (-PB*AQ+PYBIJ) / R Y + HS)*DDT 

3002 CONTINUE 
3001 CONTINUE 

4002 CONTINUE 
3000 CONTINUE 

UUIJ = 0.0 
VVV=0 .0 

IF(MAT.BQ«0) GO TO 3010 
IF (MAT.EQ. { (HAT/2) *2) ) GO TO 3012 
IF (MAT. EQ. 1) GO TO 3013 
C MAT- 3, 5, 7 

3 01 1 READ (3) iJ, V 
REWIND 3 
DO 3106 

XJ1 = X0A ( J) /2 . 0 
HT=1 . 0/(1 .0+XJ1**2) 

Hs=VHSM1 /7.0 



a i i t j i% a 1 1 i y *j f 

I>J=LV{J) 

^ DO 400S H1=1,LJ 

v - 11= LA (Ml , J) 

; IG=I. B (Ml , 3) 

>\ DO 3107 1=11, IG 

• DO 31 OR K = 1 , KM 

VL=XMJN*7 (K) 

AT J = U (I , J , K) 

BIJ= V (I , J , K) 

AU" = A IJ + TEN (I, J,K) + HS*BTJ 
AVZ=BIJ+3SN(I r‘ 7 , K) -HS*AIJ 
P1= (ATIZ+XJ1 *A VZ) *HT 
P2= (A VS-XJ1* AUZ) *ITT 
IJU(I, J, K) =P 1 
VV (I , J , K ) =P2 
UUU=T1I10+P1*P1*VL 
V V V= V7V + P2*P2*VL 
3103 CONTI DUE 
3107 CONTI NO? 

4008 CONTINUE 
3106 CONTINUE 
10 DO 8 K=1,KPt 
DO 12 

UO ( 2 , «1 , K } =00 (IE, J,FC) 

VV(2,J,K) =VV ( I M , J , K ) 

UU (THI1 , J , K) =00 (2 , J,K) 

VV (IMP! , 3 , K) =VV {3 ,J,K) 

12 CONTINUE 
8 CONTINUE 
DEAD (2) T 
REWIND 2 

IF (HAT.EQ.7) GO TO 3120 

WRITE (3) UU,VV 

REWIND 3 

DO 3999 K=1 , K M 

DO 39 93 J=JMK,JHm 

1)0 3 99 7 T=1 , IHP1 

n (j r j,k) =un (t,j,k) 

V(I,J,K)=VV(I,J,K) 

3997 CONTINUE 
3998 CONTINUE 
,3999 CONTINUE 
GO TO 1 1 
C MAT= 1 

3013 DO 3 100 

LJ= LV (J) 

XJ1 = XJA (J) 

HT=XHT (J) 

HS= XHS ( J) 

Xil J N= XH Y (J) 

DO 4003 M 1=1, LI 
ix= la (r?l , J) 

IG=LD (Ml f 3) 

DO 3101 1=11,10 
DO 3 10? K = 1 , K M 
VL=X RON Y- (K) 

• AIU=ri(I, J,K) 

BIU=V (I, 3 , K) 

A07. = A TO +TS N (I, J,K) +BS*BIJ 

AVZ=BIJ+SEN (I ,J,K) -HS*AIJ 

AI J = (AUZ+XJ1 * AV7) *HT 

BIJ= (AVZ-XJ1 * AUZ) yht 

t?U(I,J,K)=AIJ 

vv(i f j,k) =nu 

UU 0= nuu + AT J* A T L 



V V V - V V V 1) 1 J !J> XiJ V j, 

3102 CONTINUE 
£101 CONTINUE 

4 0 0 3 CONTINUE ' 

<33 00 CONTINUE 
GO TO 10 
MAT =2 ,4 ,G 

2 DO 3 103 JHM1 

X J 1 = X J A (J) /2 . 0 
HT = 1 . 0/ (1 . 0 + XJl*XJ1 } 

H S= X B S { J ) /2.0 
L.1=LV (J) 
xfiJN^xnx (j) 

DO 4 0C4 M 1 - 1 , L ,7 
1 1 = L A (HI ,3) 

IG = LB (Ml ,3) 

DO 3104 1-11,10 
DO 3105 K = 1 , K H 
VL=X NON*?, (K) 

AIJ=U (I, J, K) 

B I J= V (T, J,K) 

AUZ = A II +TE N (I , J ,K) + ItS*BIJ 
A V7 - B II + 3E N (I ,1, K) -HS*AIJ 
P1 = {AU3+ XII *AV7.) * HT 
P2= (AVZ-XJ1*AU7)*HT 
0(1, J,K) -21 
DU ( I r 0 , K) - PI 

V (T , 0 ,K) = P2 

V V {I , J , X } -P2 
tHlTJ-Umi+ n 1 *VL 

V V V = V V V + V2 * P2 *VT, 

3105 CONTINUE 

3104 CONTINUE 
4004 CONTINUE 

3103 CONTINUE 

DO 3 113 K= 1,K M 
DO 3114 1= I EH, IT 
"" DO 3 115 I - 2 , 1 UP 1 
T(I, J,K) = TT (T ,I,X) 

3115 CONTINUE 
3114 CONTINUE 
3113 CONTINUE 

GO TO 3120 
C M A T=0 

3010 DO 3 100 J=JMH,JF3fn 
XJ1 = X JA ( J) 

H T= X H T ( J ) 

HS-XHS (J) 

XHJN=XMY (I) 

LI=L V (I) 

DO 4 006 t?1 = 1,LJ 
11= LA (ill ,1) 

IG=IB (Ml ,1) 

DO 3110 I = 1 1 , IG 
DO 3111 K=1 f KK 
VL=Z (K) *XMJN 
AIJ-U (I, J,K) 

• BII- V (I , I , K) 

U =un (I, J,K) 

V(I ,3,K) = VV{I,J,K) 

ADZ=ATJ +TEN(T, J, K) +HS*BIJ 
A V 7 = B IJ + 5E N (1,1, K) -HS*AIJ 
P 1= (A UZ + X J 1*S 77) *HT 
P2= ( AVZ~XI1*AUZ) *HT 
UU (I ,0,K)=P 1 
VV (I, ,l r K) =P2 


76 


4 


VI 


II U U — II IJ11+ L» 1 ■•■V I vtf 
VVV=VVV+i>2*P2*VL 
3.111 CONTINUE 
*•’3110 CONTINUE 
>'4 £> 06 CONTINUE 
<3309 CONTINUE 
20 CONTINUE 

DO 40C7 K- 1 , K M 
DO 4008 J = JMH,,mM1 
U (2, J,K) -U { I M , 0 , K ) 

V (2, J,K) = V(T.M,J,K) 

00 (2 , J,K) =U17 (IM, J,K) 

VV (2 , J,K) = VV (IK,J,K) 

U (TMP1, J, K) =U ( 3 , J , K ) 

V (TMP1,J,K)=V ( 3 , J , K ) 

OH (T Mil ,J, K) — 13 IT (3 ,J,K) 

YV (IMP1 , J, K) = VV (3 , J,K) 

4008 CONTINUE 
4007 CONTINUE 
11 UUU=UUU*BRH 
WV=VVV*RPO 
uv“unu+vvv 

R ET U F M 
END 

SUBROUTINE CONTI 

COMMON T (147,3 9,5) , TT (147 ,39,5) ,SP(147,39) ,A? (5) ,00 (147,39, 5) , 

1VV (1 47,39,5) ,0(147,39,5) , V ( 147, 39,5) , W( 147, 39,5) ,»B( 147,3 9) , 

2VB( 1 47,3 9) ,WE ( 147,39) ,SfW( 147,39) ,0(39) ,XJR(38) , XHT ( 39) ,XHS (39) , 

3 XU Y (38) , 7. (5) ,X{39) ,TH (39) , Til ( 39) ,FJU (39) ,77, (6) , H , A 7 , A H , DT , A L, 

UG A K , Y , Y Y , R , RE , HK, 02 , .A 7, H R , 22 , AMY , A HP P , A U 2 , B3 , T3 P2 , A AA , RBB , DDT , A 1 , 
5UUU, VVV , RP H, K0, GAM2 , AL2 ,3 UN, RY, AM , AMR R , GANG AM, A I A L , BA N, KM PI , 

6MPP (1 47, 3 9) , M QQ ( 1 47, 3 9) ,IH, IN Ml ,JKH1 ,KMM1 , TKF1 , .7 M P 1 , M AT , 

7 UK, NNNN ,1MB , IMU, MATS , L (3 9) , KA ( 10,39) , KB (10,39) , LV (39) , NFND , 

8LA (10,3 9) , ID (10,3 9) , U N , T F (14 7,46) ,UN( 14 7,39) , TP ( 14 7,46) ,PN (147, 39) 
* , Cl. ( 147, 39) 

DIMENSION WWW (10) 


UV 


5 


3 = 


DO 10C0 
XJ=X (J) 

»ia 1 = 0-1 

X f? J~ X M ( J) 

XK1 = XH( JM 1) 
ST=1 . 0/ (XJ*Y) 
B=XM J*WT 
CCC-XM1* NT 


in 


LJ=L (J) 

DO 1 M1=1,IJ 
11= K A (Ml, J) 

IG=K B (M 1 , J) 

DO 1001 1=1.1,10 
M P= M P P (1,0) 

IM 1=1-1 


TF(SiP.EQ.I) GO TO 1200 


1 999 IP (MP.EQ. 11) 

GO 

TO 

101 1 

IF (MP.EQ. 12) 

GO 

TO 

1012 

IF ( M P.FQ. 13) 

GO 

TO 

1013 

IF (MP.EQ. 1 4) 

GO 

TO 

1014 

IF (MP.EQ. 21) 

GO 

TO 

1021 

IF (MP.EQ. 22) 

GO 

TO 

1022 

TF(I1P.EQ. 23) 

GO 

TO 

1023 

IF (MP.EQ. 24) 

GO 

TO 

1024 

IF (MP.EQ. 31) 

GO 

TO 

1031 

IF (MP.EQ. 32) 

GO 

TO 

1032 

IF { MP.EQ. 34) 

GO 

TO 

1034 


1033 DO 1113 K=1 r KHM1 

WWW (K)= ( (Oil (T,J,K) +im (T,J*1 ,K) -UU (I Ml ,J,K) ) /YJ 
1+ fVV ft, .1- K) + V Y (TM 1 ,,7 r K) 1 *P-VV (J , JM1 , Y ) *CCC) *p AN*7 m 



I \ 1 J LUNT1 NU h 

GO TO 199 8 


f 1<200 DO 1003 K=1,KMM1 

R * WWW ( K ) = ( ( (HU {1,0, K) + UU {I, JM1 , K) -UD (T FI 1,0, K) ”UTI (TM 1 , «TMl , K) ) /XJ 
^ [ 1 + (VV (I, J , K) 4 VV (I Ml , J, K) ) *0- (VV {T, JK1 , K) + W (IK 1, JM 1, K) ) *CCC) /2 


2)*7(K) 

|K)3 CONTINUE 
W GO TO 1993 
1011 DO 1004 K ~ 1 , K M M 1 
WWW ( K) - {- (00(1111 


, J,K) + UU (I Ml , J.M 1, K) ) / XJ+VV F) *D 


.0 


1-VV (101 , a Ml , K) *CCC) *2 (K) 
1004 CONTINUE 

GO TO 1998 


1012 DO 10 05 K=1,KMM1 

www ( K)- ( (un (i ,jmi , K> -ud (I hi , jmi ,k) ) /xj~ 

1 (VV ( I , JH1 , K) + VV (I Ml , JMI , K) ) *CCC) *2 (K) 
1005 CONTINUE 

GO TO 1 998 


1013 DO 1006 K = 1 , KM M 1 

www (K)= { (nu (X , j,k) + uu (i , am , k) } /xj+vv (i, j,k) *b-vy (i , jmi ,k) *ccc) - 

17 { K ) 

1006 CONTINUE 

GO TO 1998 

1014 DO 10 07 K=1,KM1 

WWW (K) = { (UU (X ,J,K) -UD (I Ml , J, K) ) /3CJ+ (VV (I,J,K) +VV ( I H 1 , J r *) ) *8) *7 (K) 

1007 CONTINUE 


GO TO 1998 

1021 DO 1 008 K = 1 , KMM1 

WWW (K}= ( — U U (I Ml , JMI ,.K) /7J-VV (TK1,JK1, K) *CCC) *2. 0*7 f K) 

1008 CONTINUE 
GO TO 1 998 

1022 DO 1009 K=1,KKM1 

WWW ( K) = (UU ( T , J M 1 , K) /XJ -VV (I, JMI ,K) *CCC) *2.0*7- (F ) 

1009 CONTINUE 
GO TO 1993 

1023 DO 1109 K= 1 , K MM 1 

WWW (K) — 2.0*7 ( F) * ( UU (I , J , K ) /X J + VV (I , J , K) *B) 

1109 CONTINUE 

GO TO 1998 


1024 DO 1110 K = 1 , K M M 1 

WWW (K)=2.0*7{K) * ( -IJ U ( I M 1 , J , K ) /X J + V V (IM1,J,K) *B) 

1110 CONTINUE 
GO TO 1993 

10 31 DO 1111 K= 1 , K KM1 

WWW (K)= ( (III) (I , JMI ,K) -UU (I Ml, J,K)-UTJ(IM1 ,JM1 , K) ) /XJ+VV (7PI1,J,K)*B 
1- (VV (I,jri1,K) + V V ( I M 1 , J F 1 , K) ) *CCC) *BAN*7 (K) 

1111 CONTINUE 
GO TO 1990 


1032 DO 1112 K-=1 # KMH1 

WWW (K) = ( (UtJ (I , J,K) + UU (T , JMI , K) - UU (TUI , J M 1 , K) ) /XJ 
1 - (VV (I Ml , JM1 , K) 4 VV (I, JMI, FT) ) *CCC+VY (I , J , K ) * B } *BAN*7 (K ) 

1112 CONTINUE 

GO TO 1998 

1034 DO 1114 K=1,KMM1 

WWW (K) = ( (UU (I,J,K) -00 (I Ml , JMI , K) - DTI (I Ml ,J,K) ) /XJ 
1+ (VV ( T , .1 , K ) + V V { I M 1 , J , K ) ) *3-VV ,K) *CCC) *BAN*7 (K) 

1114 CONTINUE 

• 98 WT = 0 . 0 

} DO 1 997 K = 1 , K M M 1 
W1= WT+WWW (K) 

W(I, J,K+ 1) "W 1 

wt= vn 

1997 CONTINUE 
1001 CONTINUE 
1 CONTINUE 
1000 CONTINUE 



J -> uf .* 


It ^ i U i\ .M 


END 

SUBROUTINE THEE BE 

» fc . COMMON 7(147,39,5) , T? ( 1 47 , 39 ,5) ,SP( 147, 397 , AT (5) , UU (147,39,5) , 

^ * 1 VV (1 4 7,39,5) , U (147, 39,5) , V (147, 39,5) , W ( 147,39,5) ,UF ( 147,39) , 

2VB ( 1 47, 3 9) , WE ( 147,39) , UN ( 147,3 9) ,Q(39) , X JA (3 5)., XTTT ( 38) , XHS (38) , 

3 XX Y (39) , 7 : (5) ,X(3 9) , XM (3 9) , FIT (39) ,FJU (39) ,7-5 ( 5 ) ,U, A7, A H , DT , AT,, 

1 4 G A M , Y , Y Y , R , R P , H B , H2 , A 7- E H , R2 , AHY , A HP E, A H 2 , B-B , B B 2 , A A A , BBE , DDT , A 1 , UV, 

5 IT on , VVV, RR Hr AQ, GA M2 , A 1.2 , B H U , R Y , AM , AMR F , GAMG A M, A L AT., EA N, fCf? Pi , 

6K.PP (1 47,3 9) , K QQ ( 1 47,39) , J M , JR , K M,I M M 1 , JMH1 , KHH1 , IMP1 , JKP1 , M AT, 

7 NK, NNNN ,IMB, JMH, NATS, I, (3 9 ) ,KA(10,39) ,KB(10,39) ,I,V (39) , N END , 

8LA (10,39) , LB (10,39) ,M[«,7F. (147,45) ,TJN( 147, 39) , TD (147,45) , RN ( 147, 39) 
* , CL { 147, 3 9) 

CON M CN/T EM P/TFN (14 7,3 9,5) 

RG=H *490 . 0 
IP(KAT.BE.H) GO TO 6 
NAT - 0 


73 ,' 


WRITE (6,2 34) NN 

234 FORMAT (1U ,3 H END , 4 X ,5 HTIME=I6) 
6 CONTINUE 

SBC = 8 , 3 E- 1 T* 0 . 0 1 6 6 6 


DO 2 0 J = .7 M H , J M 
JH1=J-1 


Q (I) 
X J= X ( J) 
091=1 + 1 


AA= 1 * 0/ (R.2*XJ) 
A.A2 = AA/2* 0 
XJ1 = X (INI ) 

XMJ= XF (J) 

X N 1 1 ■= X N (INI) 

BB2 = X«J*A A2/Y 
AH22=AH2* XKJ/XI 
EB1=BB2* 2.0 
CC 2 = XN7 1+AA2/Y 
CC=CC2*2. 0 


A HN “ A H Y /XI 
A HI =AHPR/ (XJ*XJ) 

AH3=AH2*XMJ1/X J 
I. J=L (J) 

DO 4 011 N 1=1, II 
II=K A (Ml, J) 

IG= K B ( M 1 , I) 

DO 21 I=II,IG 
UN=UN (I, J) *51 - 48 
TEI J-TE (I,J) 

CIN=CL (I, I) *0.125 

m p = x r p (i,i) 

RSN=QI* (1.-0. 7*C.LN) 

2 ♦*#^^CONPUTATICN OF PRESSURE*********************** 
P 7 = 0 - 0 
PZ 1=0.0 

DO 4018 K = 1 , K MM1 
KP 1 = K + 1 

P7=P7-77 ( K P 1 ) *AL2 * (TT ( T. , J , K) +TT (I, J,K Pi) ) 
PZ1=PZ1 + P 2 * 7. (KP1 } 

4018 CONTINUE 

SP(I ,0) =P71*HG 
IP1 =T+1 
IN 1 = 1-1 
UBT I=DB (T , I) 

UBM 1 J = UB (TNI ,7) 

TJBIM 1 =T.T B ( T , J N 1 ) 

VET I = VB (T , J) 

YBM 1I = VB (INI , I) 

VBIM 1 =VB ( T , J N 1 ) 



UttH III I' — I ! H { .1 I , J n I ) 

VBMim = VB (I Ml , JM1 ) 

Tft = TE (I, J) *273.15 
r ' TW= T D (I , .7) + 27 3.15 

^ j23 TS=T{I/t1 ,1)+2 73 .15 
PC B = 101 3.25*0.1 
ESS -10. 0** (8 . 4051-2353. O/TS) 

FSD = 1 0 . 0** {8. 4051-2353. 0/TW) 

PLN=SBC*TS**4*0.9B5*(0. 39-0.05*SQRT (BSD* 10. ) ) * { 1 . - 0 . 6*CLN** 2) 

BNS= RSN- RLN 

BN ( 1 , 0) “ ^ NS 

QS=0 .62 2*ESS/ ( n CB— ESS) 

QA=0. 622*FISD/ (PCB-ESD) 

DO 22 K” 1 , KM 
KP1 = K+ 1 
2Z1=Z 7, (K P 1 ) 

ZZ K=ZZ (K) 

ZK=Z ( K) 

C VERTICAL DIFFUSION (HT, HS) AND VERTICA L ADVECTIOM ( WT , WS) 

WIJ = W (T ,J,K) 

CC1 =1 .0 / (R2*ZK) 

KM1-K-1 

TTP I! =TT (IP1 , J , K) 

TTM1 J=TT(TJH1 , J,K) 

TTIP 1 =TT ( T r J? 1 ,Y) 

TTIE 1-TT (I ,J!1 1 ,K) 

TT.I J = TT ( I r I, K ) 

TIJ-T (I ,J,K) 

IF{K .FQ. 1 .AND .{ (TIJ-TETJ) +0,509 *TS * ( QS-QA) ) . GE . 0 . 0) GOTO 333 

IF(K . EQ. 1) GO TO 3 34 

TXI J = (T (I ,0, KH1 ) “II J) /77.X 

WT= (TTIJ + TT {I, J, KH1) )*WIJ 

IF(K.NE.KM) GO TO 2003 

TYI J = C. 0 

PK=0.0 

PK 1 =0.0 

GO TO 2004 

333 CH=0.033 
CO=1 .2-3*0.001 
GO TO 133 

334 CU = 0.013 
CO=C. 41*0. 001 

...133 QH=1 . 2* 0 . 001* 0. 34 *CO* U M* (TIJ-TETJ) 

QE= 1 . 2*0. 001* CO* 5 95.* BE * (QS-QA) 

TXT .1 = (PN5-QH-QE) *H/AZ 

WT= 0.0 • 

2003 TYIJ = (TIJ-T (I ,J,KP1) ) /ZZ1 
W.I J P- W (T , J , K ? 1 ) 

PK= (TTIJ+TT (T ,J,KP1 ) ) *WTJP 

2004 H S= A ZHH/ZK 

HT= (UT-PK) *CC 1 
HT=HS* (TXIJ-TYIJ) 

C HORIZONTAL DIFFUSION (HT,US) AND HORIZONTAL ADVECTION (WT,WS) 

TXIJ= (T (IP1, J, K) -TIJ) *AH1 
TXM1 J= {TIJ-T (TMl , J , K) ) * AH 1 
T YI J- (T (I, JP1 , K) -TIJ) * A H 2 2 
T YI tl 1 = ( TIJ-T ( I,JF1,K) ) * AH 3 
VIJ=VV (I , J , K) +VBI J 
VM1 J = VV (INI , J ,K) + VBM1 J 
VIM 1 = VV (I , J Ml , K ) + V B I M 1 
VM1K 1 = VV (I M 1, JMl r K) +VBM1M1 
UI J" II 0 {I r J , K) +UB.IJ 
UH1 J = DU (INI , J r K) +UBM1 J 
HIM 1 = 011 ( I, JM1 r K) MINIMI 
U Ml K1-UTJ (I Ml , J .M1, K) +UBM1M1 
TFfMP.EO. 1) GO TO 2 001 


if lin I Ij >1' ! ... '■J u o 

IF{tfP.HQ. 12) GO TO 2012 

IF ( P1P.FO . 13) GO TO 2013 7$ 

7 IF (KP.EQ. 14) GO TO 2014 
IF ( MP*FQ .21) GO TO 2021 
IF (MP.EQ. 23} GO TO 2023 
IF ( KP.EO. 24) GO TO 2024 
IF (MP.EQ. .3 1) GO TO 2031 
IF (MP.EQ- 32) GO TO 20.32 
IF (K P.FQ. 34) GO TO 2034 
IF (MP. HQ. 22) GO TC 2022 
2033 F=UIJ+ni Ml 

G-( VIJ+VU1J) * B3 2 
FF- VI fl 1 *CC2 

WT= VT + ( { (TTI J+TTP 1 J) * F - ( TT IJ + TTM 1 J) *T1M 1 J) * A A 2 + 

1 (TTU+TTI PI ) *G - (TTTJ+TTI Ml) *FF ) *BIJW 

HT= HT + (TXIJ+TYIJ- (TXM1 J+TX T Ml ) /2 . 0) *B UK 
GO TO 2005 
2001 F=UTJ+0IM1 

G= (VIJ+VM1 J) *BB2 

FF = I3M1J+UB1*n 

GG- ( VIH1 + VM 1M 1) *CC2 

WT= VT + ( (TTI J+TT PI J) *F - (TTIJ+ TTM1J)*F? ) *AA2 

1+ (TTIJ+TTIP1 ) *G - {TTI J + TT IM 1) *GG 

.HT- T XI J-T X M 1 J +TY I J-TYIM 1 4 HT 
GO TO 2005 

2011 F=(UM1M1+ IJK1 J) *AR 
G=V.M1 J*BB 1 
FF=VM1H 1*C C 

WT- WT- (TTTJ+TTM1 J) *F + (TTIJ + TTIP1) *G 

1- (TTTJ + TTI Ml ) *FF 
HT=- TXM 1.1*2 *0 +TYT J-TYIM1+ ITT 
GO TO 2005 

2012 F=(VIMl+VMl Ml ) *CC 

$T= $T + ( (TTTJ+TTP1 J) *HIM 1“ (TTIJ+TT Ml J) *TIM1M1) * A A 
1- (TTIJ+TTTM 1) * F 
HT=T XI J-TXF1 J-TYIM 1 *2 .0 + HT 
GO TO 2005 

2013 F= (0TJ+UTM1) * A A 
G -VI J*BB 1 

FF— VI Ml * CC 

WT= WT+ (TTIJ +TTP1 J) *F + (TTI J+TT TP 1) *G 

1- (TTIJ+TTIH1) *FF 
HT=TXIJ*2 .04 TYI J-TYIM 1+HT 
GO TO 2005 

2014 F= (VIJ+VMU) * BB1 

WT=WT+ ( (TTIJ + TTP1 J) *0IJ- (TTIJ+TTMU) *0M1 J) *AA 
1+ (TTI J+ TTI Pi ) *f 
HT=TXIJ-TXM1J+TYIJ*2- 0+HT 
GO TO 2005 

2021 F=TJ M 1 Ml * A A 
G = VM 1M1*CC 

WT- WT- (TTIJ+TTM1 J) *F * 2 .0- (TTI J 4-TTIM1) *G *2-0 

HT= HT- (TXM1J+TYTM1) *2.0 
GO TO 2005 

2022 ?=0IM1*AA 

G=VIM1*CC 

»T=WT+ ( (TTI J+ TTP1 J) *F -(TTTJ+TTJM1) *G )*2.0 

HT= (TXI J-TYIM 1 ) *2 « 0 + H T 
GO TO 2005 

2023 F = U I J * A A 
G=V I J*EB1 

WT= WT+ ( {TTIJ + TTP1 J) *F + (TTI J + TT IP 1) *G }*2.P 

HT= (TXI J + TYJJ) * 2 , 0 + HT 
GO TO 2005 
202 4 F-UM 1 J* A A 



) * 2.0 


WT= W T“ ( (TTIJ+TT Ml J) *F - (TTIJ+TTIPI) *G 

HT= (TYIJ-TXM1 J) *2 . O+HT 
GO TO 2005 
►2031 F=UMl2i1 + UMl J 
G=VH 1J*BD2 
PF= (VIP, H VM1M1) * CC2 

WT= WT+ ( ( (TTJ J+TTP1 J) *01 Ml - ( TTIJ+TT M 1 J ) * F 
1 (TTIJ +TTIP1) *G - (TTIJ + TTIMI) *FF 

HT= ( (TXI J + TY.T. 0) /2 ,0-TXM IJ-TYIMI ) *BUN+ HT 
GO TO 2005 
2032 F = (! I M 1 + u I J 
G=V I J*BB 2 

FF= (VIM1+ VMlMl) *CC2 
WT= WT+ ( ( (TTIJ+TT P 1 J) * ? - (TT'IJ +TTM1 J) 

1 (TTIJ+TTIPI ) *G - (TTI J+TTIM 1) *FF ) 

HT= ( (TYIJ-TXM 1J) /2.0+TXTJ-TYIM1) *BUN+ HT 
GO TO 2005 

2034 G- (VIJ+VH1J) * BB2 
F=VM 1 M 1* CC2 
FP-UBlJ+Uin!11 

ffT=MT+( ( (TTTJ+TTP1 J)*IJIJ- (TTIJ+TTH1J) ) 

1 (TTIJ+TTIPI) *G - (TTIJ + TTIMI) *F ) 

HT= ( (TXI J-TYI FI) /2.0-TXM1 J+TYIJ) *BUN+ HT 
2005 TEH { I , J , K) = D P T * (HT-HT/O . 93) 

22 CONTINUE 
21 CONTINUE 
4011 CONTINUE 
20 CONTINUE 

IF(MAT.EQ.Q) GO TO 20 36 
IF (KAT.EQ. ( (MAT/2)*2) ) GO TO 2037 
IF {MAT. EQ- 1) GO TO 20 3 8 
: MAT=3,5,7 

203 5 HEAD (2) TT 
REWIND 2 
DO 3 8 J=JKH,JM 
LJ^L (J) 

•" DO 4 014 fi1=1 f LJ 
II-KA {Ml , J) 

IG= KB (M 1 , J) 

DO 39 I -I I , TG 
DO 40 K= 1, KM 

AT (K) ~TT (I, J, K) +TEN (I , J,K) 

40 CONTINUE 
CALL ADJUST 
DO 41 K=1,KM 

TT (X , J, K) =AT(K) 

41 CONTINUE 
39 CONTINUE 

4014 CONTINUE 
38 CONTINUE 

IF (MAT. EQ.7) GO TO 2500 
9 DO 2 K~1 /KM 
DO 3 J = J M H , J M 
TT(TMF1 # J,K)=TT(3,J,K) 

TT (2 , J, K) ~TT (I M/ J, K) 

T ( 2 , J / K ) - T (IM/J/K) 

T (IMP1/J/K) =T (3/J/K) 

3 CONTINUE 
2 CONTINUE 
WRITE (?) TT 
REWIND 2 
GO TO 10 
; MA T= 1 

2038 DO 30 J-JMH/JM 
I.J-L (0} 


) *AA? 
) * B IT N 


* A A 7 . + 
*BDN 


* A A 2 + 
*BUH 





T 



DU 4 U 1 2 W t - I , LJ 
II=KA(M1,J} 

IG= KB (HI , J) 

DO 31 I— I I , I G 

DO 32 K=1,KK 

AT(K) -T (I # ,J,K) + TE N (I, 

CONTI NOE 

CALL ADJUST 

DO 33 K=1,KH 

T(I,J,K)=TT(T,J,K) 

TT (I , J , K) -AT ( K) 


33 CONTINUE 
31 CONTINUE 
4012 CONTINUE 
30 CONTINUE 
GO TO 9 


J, 




C HAT-2, 4,6 

2037 DO 34 

LJ-L (J) 

DO 4013 M 1 =1 , L J 
11“ K A (Ml , J) 

IG-KB (K 1 , J) 

DO 35 1-11,10 
DO 3 6 *0=1, KM. 

AT(K)=T(T,J,X) +TEN (I, J, K) 

36 CONTINUE 
CALL ADJUST 
DO 37 K-1,KH 
TT(I, J,K)-A?(K) 

37 CONTINUE 
35 CONTINUE 

4013 CONTINUE 
34 CONTINUE 
GO TO 2500 
C MAT-0 

2036 DO 42 J-JHH, JM 
LJ-L (J) 

DO 4015 HI-1 , L J 
IT— KA (Ml , J) 

IG-KB (Ml , J) 

DO 4 3 I — II, T G 
DO 44 K-1 , KM 

AT ( K ) -T ( I , J , K ) +TE N (1, J, K) 

44 CONTINUE 
CALL ADJUST 
DO 45 K-1, KM 
T(I,J,K)-T?(T,J,K) 

TT (I,J, K) — AT (K) 

45 CONTINUE 
43 CONTINUE 

4015 CONTINUE 
42 CONTINUE 

2500 CONTINUE 

DO 4016 K-1, KM 
DO 4017 J - JMH , JM 
T{2 , J,K) -T (IM ,J,K) 
TT(2,0,K)-TT ( I M , J , K ) 

T (IMP 1 , J , K) =T ( 3 , J , K ) 

T T ( I M F 1 , J , K ) - T T ( 3 , J , K ) 
4017 CONTINUE 

4016 CONTINUE 

10 DO 4019 J-JMH , JM 
SP(IMF1,J)-SP ( 3 , J ) 

4019 CONTINUE 
RETUB N 


END 



A 1*1 

T) b 


bUh HU Ui 1 N Hi 

RETURN 

END 

SUBROUTINE 

RETURN 

END 

SUBROUTINE 

RETURN 

END 

SUBROUTINE 


Eh l.I' 


GFOST 


PREP A 


SHADE { ICTL, RCTL, BVAL) 


. LT N { 1 3 2 ) , RUN (132) 


REA I RCTL (2), BVAL (44), 

INTEGER ICTL (8) 

LOGICAL* 1 BLANK/* '/, LP/'PA RP/M 1 / 

LOGICAL*! FRSTLN, RATAL, NS HADE, CORD, CLIN ( 137) 


CTRL (2 0) 


£***#****$*£ #*****$**«*:^< *£******#£$ ******* ******** ********* ********<< 
C CALL SHADE (ICTL, RCTL, R FI ELD) 

£4****$** $********$ ************************ ******** ******* &********** 


c 

ICTL ( 1) 

4 • « 

CALL CODE 

c 


= 0 

-NORMAL INTRA-FIELD CALL 

c 



(ONLY F FIELD IS USED, ICTL (2) THRU ICTL (8) AND 

c 



RCTL APE IGNORED) 

c 


= 1 

-BEGTNNING-OF-FXELD CALL 

c 



(ONLY ICTL AND RCTL PARAMETERS ARE USED, F FIELD 

c 



IS IGNORED) 

c 

ICTL (2) 


OUTPUT TAPE NUMBER 

c 

ICTL (3) 

♦ * » 

LEFT (SHADE) MARGIN 

c 

ICTL (4) 


NUMBER OF SHADING LINES 

c 

ICTL (5) 

* « -a 

NUMBER OF FIELDS 

c 

ICTL (6) 

m m m 

FIELD WIDTH 

c 

ICTL (7) 

• a a 

DOES THE CALLING PROGRAM PRINT INTERSPERSED DATA 

c 


= 1 

-NO 

c 


= 2 

-YES 

c 

r 

ICTL (9) 

■* * * 

IGNORED 

C 

RCTL (1) 

«*• * m 

CYCLE POT Mr 

C 

r* 

" RCTL (2) 

« » a 

SHADE INTERVAL 

C 

R FI ELD 

« • « 

ONE DIMENSIONAL VECTOR WITH NUMERIC INFORMATION 

C 



TO BE SHADED 


£********** ** * * * * * * * * * * **** ** ************** ** * * * * * * ** * * * * $ 
IF (ICTL ( 1 ) - NF. 0) GO TO 1000 

X L~ {BVAL ( 1) -CP) /SINT . 

M2=LM 

IF (CCND) GO TO 2100 
. BLIN (LM) = {XL - At. JM (IH) ) /VDIVS' 

3021 DO 3029 K=L'F , NF, I SKP 
M1=M2+1 


** ******* if. * 


H2=M2+IFS 

XP = ( (BVAL (K) - CP) /SINT- XL) /FW 

DO 3029 J=M1,M2 

XJ.-XL+XR 


3029 BLIN (J} = (XL-ALIN (J) ) /VDIVS 
3051 DO 3059 L=1,NSL 
3041 DO 3049 K=LM,HR 

A IT N ( K) = A LI N ( K) +BLI N (K ) 

TBLC=AHOD (ALIN (K) , TIE N) +1 , 

IF (TBLC .IT. 1.) TBLC = TBLC+TLFN 
it-tbic 

3049 CLIN (K+2) =CTIU. (IT) 

3 059 WRITE ( NT P, CL IN) 

IF ( * NOT. DATA L) RETURN 
3061 DO 3069 K*LH, NR 
3069 A LI N (K) = ALIN ( K) +BLI N ( K) 

R ET U R N 



00000100 
00000 120 
00 0 001 4 0 
00000160 
00000190 
00000200 
0000022 0 
000 00 240 
00000260 
00000290 
00000 300 
00 0 00 3 20 
00 0 00 34 0 
0 00 003 6 0 
000 00 390 
00000400 
00 0 00420 
000 0044 o 
00000460 
0 0 0 00490 
0 0 0 00 5 0 0 
00000*2 0 
00 0 00 c 4 0 
00000*60 
000 00*9 0 
00000 600 
00000620 
000 00 64 0 
000 0066 0 
000 006 bo 

00000700 
00000720 
000 0074 0 
0000076 0 
00000780 
OOOOOflOO 
000 0087 o 

000 00 54 0 
000 00 86 0 

000 00830 
00000°00 
00000920 
00000940 
00000960 
00000980 
00001 000 
00001020 
00001 040 
0 00 01 06 0 

00001 080 
00001 100 
00001 120 

000 01 14 0 
00001160 
00 0 01 18 0 

00001 200 
000 01 22 0 



K> - U VJ 1 1 


Jiuu xr (N SHADE) uu i:o 2 uuu 
COND=. FALSE. 

ALIN(1M)=XL 

*2121 DO 3129 K=LP,NF,ISKP 
"”7 M 1 = M 2+ 1 

H2=M2+JFW 

XR = ( { EVA L (K) -CP) /SIWT-XL) /FW 
DO 3129 J = 

XL— XL+XR 
312 9 A LI N { J) = XL 

IF (RATAL) RETURN 
3141 DO 3149 X=LM,NR 

TBLC= AMOD (AL IN (K.) , TIEN) +1 . 

IF (T13LC .LT. 1.) TBLC-TB LC+T I.EN 
I T= T B LC 

3149 CLIN ( K + 2) =CTBL ( IT) 

WRITE ( N TP, CL IN) 

RETURN 

2000 IF (FRS7LN) GO TO 2030 
2011 DO 2019 L = 1 , NSL 
2019 WRITE (NTT?, CLIN) 

2030 FFSTLN=. FALSE. 

RETURN 

1000 CONTINUE 

I5KP-ICTT. (1) 

LF = I SKP+ 1 
NTF = ICTL{2) 

LM=ICTL (3) 

NSL=IC?T,(4) 

NF-ICTL (5) 

IFW=ICTL (6) 

DATAL=ICTL (7) * NE. 1 
VDI V S=MSL 

IF (D AT AT.) VDU/S=VDIVS+ 1. 

IF (LM.LT.2) L»*2 

IF (IFW.LT.2) IF W=2 
FW=TFW 

NS- (NF-1) /ISKP+1 
MR*T,M+ (NS-1) *IFW 
IF (NR. LE. LI) GO TO 1011 

NF= (1L-LM)/IF»*XSKP+1 

NS= (NF-1 ) /ISER+1 
KR=LM+ (NS-1) * 1FW 

101 1 CP-ECTL (1) 

SI N T= RCTL (?) 

FRSTL N= . TRT1E. 

CON D=. TRUE. 

NSHADE=SINT.EQ.O. 

IF {. NOT. NS HADE) GOTO 1021 
SIM T- 1 . 
c?~o . 

MR = 1 
1 . 11-1 

1021 DO 1029 1=1, LK 

102 9 CLIN (I f 2) = BLANK 

CLI N (1) =LP 
CLIN (2) = QM 
CLIN (MR+3) =QM 
CLIN (HR+4) =RP 
RETURN 

DATA LL/132/ 

DATA TLEN/20.0/, CTBJ./1H0,1H ,11.11,111 , 1H2,1H , 1E3,1H ,1H4,1H 
*, 1 H 5 , 1 H ,1(56,111 , 1 H 7 , 1 H ,1HB,1H ,1H9,1H / 

END 



00001 260 
0 00 01 2 0 
00001300 
00001320 
00001 340 

noo 01 36 0 

000 01 380 
00001400 
00001420 
000 01 44 0 
0 00 01 4 6 0 
000 014 80 
0 0001 ECO 

0 00 01 92 a 

00001 54 0 

000 01 560 
00001 580 
0 0 0 01 6 0 0 
00001 620 
00001 640 

00001 66 0 
16 80 

00001900 
000 01 77 0 

0 00 01.74 0 
00001760 
00001780 

00001 800 
000 01820 
000 01 840 
HOO 01 98 0 
00001 380 
00001900 
00001920 
000 01 94 0 
00 0 01 96 0 
0 00 01980 
00003000 
0 00 02 02 0 
000 02 G 4 0 
0 00 02060 
0 0 0 02 0 9 0 
00002100 
00002 120 
00002140 
000 021 6 0 
0 0 0 02 1 8 0 
00002200 
00002220 
0 0 0 02 2 4 o 
000 022 8 0 
00002239 

0 0 0 02 3 0 0 
0 00 02 32 0 
00002340 
00 0 02 36 0 
000 03 38 0 
0 0 0 02 4 0 0 
00002420 
00002440 
0 0 0 02 4 6 0 
0000/ 480 
00002500 
00002530 


nn,OT 


subroutine prep 


DATA: jmi 

OCEAN SHAPE (SUBROUTINE SHAPE?) 

H, 2(1} 

A 2, A H # A !T 

Vi E / H W (GIVEN FROM CARDS, OH A TAPE OR THE SUBROUTINE WIND- 
0 (J) (HEATING FUNCTION, REFERENCE ATMOSPHERIC TEMPERATURE) 

AT (K ) (INITIAL IE HP ERA TO HE) 

TIB, V B (GIVEN FRO M CARDS, OP. A MAGNETIC TAPE, OP. THE SUBPOUTIN' 
BA UO) 

BB, A AA, R BB (AAA: GRID SIZE IN LATITUDE, BBB : GP ID SIZE IN 
LONGITUDE. BOTH IN DEGREES. BB IS THE LATITUDE CORRESPOND- 
ING TO ,7»0.) 

COMMON r (H7 f 19,5) ,TT (1 47,39,5) ,SP{14 7,39) , AT (5) ,00 (147, 39, 5) , 

■m ( 14 7, 3 9,5) , 0(14 7,39,5) ,V (1 47,39,5) , W (147,3 9,5) ,UB (147 ,39) , 

2VB{1 47,3 9 ) ,WF (147 ,39) , ON (147 , 39) , 0 (39) , XJA (38) , XHT (38) , XUS ( 38) , 

3 XMY (30) , 2(5) ,7. (39) , XH (39) ,FU(39) ,FJU(39) ,ZZ (6) , B , AZ , A H ,DT , A L, 

4 G AM , Y , I Y , R , R R , II H , H2 , A ZH H , R 2 , A H Y , A HR P , AH2 , B B , B B2 , A A A , BE B , DDT , A 1 , 0 V , 

5T1DU , VVV,PPH, A 0, G A M2 , AL2 , BO N , BY , AM, AMP F , G AM G A K, AI. AL, BA K , KKP1, 

SSPP (147 ,39) ,MQO (1 47,39) , I K , J M , E M, IN Ml , JMM1, KH » 1 , IMP 1 , IN? 1 , M AT , 

7 NK, NNNN, IMK, JMH, MATS, L (39) , KA (10,39) ,13(10,39) ,LV ( 39) , NEND , 

8LA ( 10,39} ,LF ( 10,3 9) 

READ (5,3 0) JMH 

80 FORMAT (515) 

JMP1-JM+ 1 
JMM1-JH-1 
XHPl=IK+ 1 
IHfil-IH- 1 
KHH1-KM-1 
CALL SRAP52 
CALL SHAPE2 
DO 4 7 -K-1 , KM 

.... DO 4 8 J=1,JKP1 
DO 4 9 1-1, IMP 1 
Tin (I, J, K) =0.0 
VY (I , J, K) =0.0 
W (I , J,K) =0.0 
U (I , J,K) =0.0 
V (I , J,K) =0 .0 
TT (I , J, K) =0. 0 
T (T , I , K) =0.0 
49 CONTINUE 
48 CONTINUE 
47 CONTINUE 

READ (5, 8 1} If, (7 (1) ,1=1, KM) 

81 FORM AT (El 5.7, 5? 10 . 2) 

SPITE (0,13) ( Z (1) ,1=1 , KM) 

13 FOP MAT ( 1 H ,19 HDEPTHS FOP. T, U, V=5F10.5) 

KMP1=KM+1 
ZZ(1) =-2. 0*3(1) 

DO 82 1=2, KM 
ZZ(T) =2 (I-D - 2 (I) 

13 2 CONTINUE 

ft 22 (KMP1) = 2.0* (1 .0+2 (KM) ) 

DO 83 1=1 , KM 

7, (I) = (22 (I) +22 (1+ 1) ) /2. 0 
83 CONTINUE 

WRITE (6 , 33) H, (2(1) ,1 = 1 , KM) 

33 FORMAT (1H , 2H H=E 1 5 . 7, 2X , 2 H 2=6E1 5. 7) 

WRITE (6,21) (22(1) ,1=1, KM PI ) 

21 FOR M AT (1 H , 3 HZZ=7 El 5. 7) 



KnAD(b,HU) h'C- , A a , A 15 
*34 F03N. AT(3E15.7) 

4 c i67 T)0 4 568 J = 1 ,39 
**V DO 4569 I=1,IHP1 
* ? WE (I, 3) =0.0 
*1- WN(I,O)=0.O 

CONTINUE 
4W8 CONTINUE 

CALL WIND 
CALL WIND 
AL= H/A7» 

DO 2 0 00 ■ J- 1 , 7 HP1 
DO 2001 I- 1 , 1 MP 1 
WE (I, J) =WE (I, J) *AL 
MP=HPP(I, J) 

IF(MP.LE.I) GO TO 200 3 
IF ( HP. (3 T. 5) GO TO 200 3 
DO 2 0 02 K=1,KM 
T (I, J,K) =0.0 
TT(I , 7, K) =0.0 
W (I , J ,K) -444. 444 

2002 CONTINUE 

2003 IF { BQQ (I , J) . ME. 2) GO TO 2001 
DO 2 0 04 K= 1 , K M 

0 (I, J,K) =3 33. 33 3 
UU (I, 0, K) =333.333 
VV(I,J,K) =444.444 
V (I, J,K) =444. 444 

2004 CONTINUE 
2001 CONTINUE 
2000 CONTINUE 

HEAD {5, 5) (0 (J) JM) 

5 FORE AT(1 2F5- 0} 

DO 5094 J-3,38 

REA R (5, 50 95} (AT (K) , K= 1 ,KH) 

509 5 FQF M AT ( 5 F 1 0 . 2 ) 

L J= L (J) 

DO 5 096 LL=1,LJ 

1 tl 1 = K A { L L , 0 } 

IP1 = KB (LL, J} 

DO 97 I=IH1,IP1 
DO 98 K = 1 , K M 
T(I, 0,K) = A?{K) 

TT (I ,0, K) = AT (K) 

98 CONTINUE 
97 CONTINUE 
5096 CONTINUE 
5094 CONTINUE 

DO 200 1=3,38 
DO 201 K= 1 , KM 
T(IKP1,J,K)=T (3, J , K) 

TT (IMP1, J,K) = TT(3,J,K) 

T(2,J,K)=T(IH,J,K) 

TT { 2 , J, K ) =TT ( IN , I , K) 

201 CONTINUE 

200 CONTINUE 

C CALL BA HO 

CALL BARO 

C AL= COEFFICIENT OF THERMAL EXPANSION 

3 AL=0.2E-3 

ALAL=AL*2.0 

C HEAT FLUX ACROSS THE SEA SURFACE IS SPECIFIED AS 70.0*(Q(J)“ 

C 11(1,0,1}) CAJ./DAY 

G AK GAH=H*7 0.0/(AZ*0.8 64E5) 

B=6.37E8 

READ (5, 1) BB, AAA, BDB 




1 FORM AT(3F10.0) 
Hn=3. 1415926/180- 0 


c > R Y= AAA* HII 

t r ' Y-HH* EBB 

* r 


5 



BO 10 I~1 r JMP1 
A 1= (BB+ BBS* FLOAT (I) } *HH 
X (I ) = RY* COS (A 1 ) 

0 CONTINUE 


CD=1 4 .58 E-5/ ( Y* R Y ) 

BO 9997 0 = 1,38 
XM (0) = (X (J) *X (J*1) ) /2 .0 
PJU ( J) = CD* (X { J) “X (J + 1 ) ) 
9997 CONTINUE 
H 2= H* 2 . 0 
YY- Y* Y 

C AAA AND BBB ABE HEIGHTS 

C OSCILLATION 

AAA =0 .6 
BBB = 0 .4 


RR=R*R 


RRH=RR*H 
AL2 - AL/2 . 0 
HH=H*H 
A2HH=A2/HH 
R2-R*2. 0 


IN I BE IMPLICIT SVUEME FOR THE INERTIA 


AHY=AH*Y 
AQ=24 5 . 0*H 
BB=YY*RR 
A 1= AM/BB 
AHK R= AH/PR 
BUN = 4 * 0/ 3.0 
BAN =2 „ 0/3 . 0 
AMR R = AM/EI? 
hH2~ AH/BB 
RY=R* Y 
BB= 1.0/ (Y*R) 

BB2 = BB/2 . 0 

RETURN 

END 

SUBROUTINE EARO 

RETURN 

END 

SUBROUTINE WIND 

RETURN 

END 

SUBROUTINE SHAPE2 

COMMON T{ 147, 39,5) ,TT( 147, 39,5) ,SP (14 7, 39), AT (5) , DO (147,39,5) , 

1VV (147,39,5) , U (14 7,39,5) ,V <147, 39,5) , W (147, 39,5) ,BB (147,39) , 

2VB{ 147, 39) , WE (147,39) , WN ( 1 47,39) ,Q (39) ,XJA(38) ,XHT(38) ,XHS (38) , 

3 XMY (38) , 2(5) ,X (39) , XM (3 9) , FO(39) ,FJU (3 9) ,Z2 (6) ,H , &2, AH,D? , AL, 

4 GAM , Y,YY ,R,RB ,HH, H2,AZHH, R2, AHY , A HER , AH2 , BB , BB2 , AAA , BBB , DDT , A 1 , UV , 

5 UUU , VVV,RRH, AQ, GA M2, AL2 ,BUN,BY, AM , AMR B ,GAHGAM, ALAL,BAN, KMP1, 

6 MPP (1 47, 3 9) , M QQ (1 47,3 9) ,IM r JM,KK,IMKl ,JMM1,KMM1 , I MP1 ,JMP1,NAT, 

7 NK r NNNN, IMH, JMH, MATS , L (3 9 ) , K A ( 1 0 , 39) ,KB ( 1 0 , 39) , LY (39) , REND , 
8LA(10,39) , LB (10,39) ,NN 

SUBROUTINE TO BE USED WHEN THE OC2AN SHAPE IS CHANGED. 


r 

( MPP 

: MPP-1 IE the FOUR NEAREST SURROUNDING POINTS ARE WATER-POINT , 

C » 2 ON THE LAND (EXCLUDING THE COAST LINE) . 

C FOR THE COAST-POINTS, REFER TO FIG.1. 

C 1 4 

C 34XXXXXXXXXXX XX X3 3 

C XXXXXLANDXXXXX 

C XXXXXXXXXXXXXX 

r XX?3 3 1 XX 



C 

t " 



c~ 4 

c 

c 

( 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


I I A A AA. 1 j 

XX XX 

XX 2 3 2 4XX 

xxxxxxxxxxxxxx 
3 1XXXXXXX XXXXXXX3 2 
12 

FIG.1 


KQQ 

MQQ = 1 IF THE FOUR NE ARES R S URE OU HDIN G POINTS ARE WATER-POINT, 
= 2 ON THE LAND, 

IN OTHER CASES, REFER TO FIG .2 * 


11 12 

xxxxxxxxxxxxxx 

XX XX XL A NDXXXXX 

xxxxxxxxxxxxxx 


XX 1 22 1 
XX 

21 1 1XX 
XX 

XXI 122 


221 1XX 
XX 

XXI 121 
XX 

211 2XX 


XXXXXXXXXXXXXX 

XXXXXXXXXXXXXX 

12 11 


FIG. 2 


E A ROC L.I NIC PART ON LX . LV {J} , LA (I, J) , AND L3(I,J) SHOULD BE DEAD 
FROM DATA CARDS. THEN, L (J) , KA (I, O) , KB (I ,0} , MQQ(J,0) AND 
N P P ( 1 , 0 } ARE DETERMINED IN THIS SUBROUTINE. 

THE DATA LV, LA AND LB USED HERE ARE IDENTICAL TO TEE DATA LV, 

LA AND LB TO BE USED FOR THE EABOTROPIC PART, BUT IT IS NOT 
THE 'CASE W T' T H L f KA AND KB. I, EA AND KB ARE DIFFERENT FROM 
EACH OTHER IN BOTH PARTS. 

THP1=I«+1 


jHPi-ja+i 
J MH - 3 
IMR=3 

DO 30 0=1,01$ 

DO 31 1=1, IMP 1 
MQQ (1,0) - 2 

31 CONTINUE 
30 CONTINUE 

DO 32 J=3 , JHM 1 
READ (5,3 3) M 

33 FORMAT (I 2) 

READ (5,34) (LA (T,J) , LB (I , J) , 1 = 1 , M) 

34 FOR MAT (7 (214, 2X)) 

DO 35 1=1, M 

1 1 = L A (1,0) 

I G=LB (I, J) 

DO 36 N=1I,IG 
A MQQ(N,J)=1 
^56 CONTINUE 

35 CONTINUE 
L V (0) =H 

MQQ(1,0)=MQQ(IMM1,J) 

M QQ (2,0) ~MQQ ( IM , 0 ) 

MQQ(IHP1,a>-HQQ(3,J) 

32 CONTINUE 

DO 1 0-2, JH 


S3 



J M I = J ” I 

DO 2 1-3 , m 
i IM1-I-1 

IF(HQQ(I, J) *EQ*1) GO TO 3 
* IF(HQQ(IH1, J) .EQ. 1) GO TO 3 

IF (MQQ (I Ml , J 231} - EQ. 1) GO TO 3 
IF (MQQ(I, JM1) .EQ. 1) GO TO 3 
HPP (I,J) ~2 
GO TO 2 

3 HPP (I,J) =1 
2 CONTINUE 

HPP (1 ,0} = HPP { IMH1 , J) 

HPP (2 ,J) = HPP ( I M , J ) 

HPP (IMP1, J) = M PP (3 , J) 

1 CONTINUE 

DO 4 J=3,JMM1 
H M=LV (J) 

WRITE (6 , 9) J, LV (J) , (LA(I,J) , LB (I , J) ,1 -1 , MM) 
9 FOR MAT ( 1H ,2015) 

4 CONTINUE 

DO 11 1=3, If* 

WRITE (6,12) I, (MPP(I, J) , J=3, JH) 

WHITE (6 , 13 ) (MQQ(I,J) , J=3,JM) 

12 FORMAT ( 1 H ,13, 2X, 3813) 

13 FORM AT ( 1 a , 6X ,381 3) 

11 CONTINUE 

DO 18 J = 3 , J M 
H M= 0 

MP=MPP (3, J) 

IF (HP.EQ.2) GO TO 70 

11=3 

MM= HH + 1 

7 0 DO 5 0 1= 4 , IKM 1 
HP PI PP (I, J) 

IF(MPFIJ.EQ. HP) GO TO 60 
GO TO ( 16 , 20 0) , MPPIJ 
16 II=I 
" MM + 1 
GO TO 60 
200 KA (MM,J) =11 
KB ( M M , J) = I- 1 
60 MP= MPPI J 
50 CONTINUE 

IF (HPP (IM, J) * EQ. 2) GO TO 80 
KA (HM,0) =11 
KB(MH,J)=IM 
00 L (J) = MM 

WRITE (6, 9) 0, L (J) , (KA (I,J) , KB (I , J) ,1=1, HM) 
18 CONTINUE 
DO 8004 

JP1=J+1 ' 

«1M1 = J-1 

DO 80G5 I=IMH,IM 

IF(MPP(I,«3) . B Q . 2 ) GO TO 8 005 

IP1 =1 + 1 

IM1=I-1 

K=0 

IF (HQQ(I,0) .EQ.2) K=K + 1G00 
IF (HQQ(TK1 ,J) .EQ. 2) K=K+100 
IF (MQQ (I Ml, Jill) - E Q. 2) K=K+10 

IF {MQQ (I, JM 1) - EQ. 2) K=K + 1 
IF (K - EO* 0) GO TO 8005 
IF (K- EQ. 11) GO TO 42 
IF (K.EQ. 1 1 00) GO TO 1100 
IF (K . EQ. 1 10) GO TO 110 
IF(K *E0* 1001) GO TO 1001 



IF (K.EQ. 1110) GO TO 1110 
' IF (K.EQ. 1 101) GO TO 1101 
V; IF (K.EQ. Ill) GO TO 111 

IF (K . EQ. 1011) GO TO 1011 
IF (K. EQ. 10) GO TO 40 
I F { K . EQ . 1) GO TO 41 
IF (K . EQ. 10 0 0) GO TO 100 0 
IF (K.EQ. 100) GO TO 100 
1 00 MPP (I, J) =32 
GO TO 8005 

1000 M PP (I,J) =31 
GO TO 8005 

41 MPP (I,J) =34 
GO TO 8005 
4 0 MPP (I,J) =3 3 
GO TO 8005 
1011 MPP (I , J) =24 
GO TO 8005 
1 1 1 MPP (I,J) =23 
GO TO 8005 
1101 MPP (I, J) =21 
GO TO 0005 
1110 MPP (1,1) =22 
GO TO 8005 

1001 MPP (1,0) =11 
GO TO 8005 

1 10 MPP (I, J) = 13 
GO TO 8005 
1100 MPP (I,J) =12 
8005 CONTINUE 
8004 CONTINUE 

DO 7501 J=2, J MM 1 
JP1 =J+1 
,7M1=J-1 

DO 7 5 CO 1=2, IM 

" ' IF (MQQ (I, J) .EQ.2) GO TO 7500 
IP1=I+1 
IM1 =1-1 
K=1 111 

IF(MQQ(IP1, J) .EQ. 2) K = 1000 + K 
IF (MQQ (I, JP1) -EQ. 2) K = K + 100 
IF{MQQ(IM 1, J) -EQ. 2) K=K + 10 
IF(MQQ(I,JMl) .EQ. 2) K=K + 1 
IF (K.EQ. 11 1 1) K= 1 
MQQ (I,J) ~K 

7500 CONTINUE 

MPP (1 , J) =MPP ( IMM1 , J) 

M QQ (I/O) - MQQ (IMM1 ,J) 

MPP (IMP 1 , J) =MPP (3 , J) 

MQQ (IMP 1,J) =M QQ ( 3 , J) 

MPP (2,0) = MPP ( I M , J ) 

MQQ (2 , J) = MQQ ( IM, J) 

7501 CONTINUE 

MPP ( 1 ,J M) = MPP (IMM 1 , JM) 

MPP (2,OM) = MPP (IM, JM) 

MPP (IMP 1, JM) = MPP (3, JM) 

DO 90C0 I=1,IMP1 
WHITE (6,12) I,(RPP{I,J) ,J=3,JM) 
WHITE (6 , 13) (MQQ(I,U) ,J=3,JM) 
9000 CONTINUE 

RETUHN f 

END 



BQ, 

l SUBROUTINE GEOST 

•» ;■ COMMON T (147,39,5) , TT ( 1 47 , 39 , 5) ,SP { 1 4 7, 39) , AT {5 ) , UU ( 1 47 , 39, 5) , 

1VV (1 4 7,39,5) ,U (147,39,5) , V (147, 39,5) , W (147,39,5) ,UB{ 147,39) , 

2VB ( 147, 39) , WE (147,39) ,WN( 147,39) ,Q{39) ,XJA (38) ,XHT(38) ,XHS{38) , 

3 XM Y (38) , Z (5) ,X(39) , XH (39) , FU (39) ,FJU(39) ,ZZ (6) ,H, AZ,AH,DT,AL, 

4GAM, Y,YY,R,RR,HH, H2 , AZHH, R2, AHY , A HER , AH2, BB, BB2 , AAA, 3BB , DDT , A 1 , UV, 

5UUU, V VV, RR H, AQ, GA M2, AL2, BU N,RY, AM, AMR R , G AMG A M, ALAI, BA N, KM PI,, 

6MPP (147,3 9) ,NQQ (1 47,3 9} ,1 H, JM,KM, IMM1 , JM M 1 , KM Ml , IMP 1 , JMP1 , ft AT , 

7 NK, NNNN/IMH, JMH, MATS,L(3 9) , KA (10,39) , K B ( 1 0, 39) , IV (39) , NEND , 

8LA{1 0,39) ,LB (10,39) ,NN 
HG=H*490- 0 
JJ=0 

DO 20 J= JMH, J M 
L J-I (J) 

DO 4011 M 1 = 1 , LJ 
11= K A (Ml , J) 

IG=KB (M 1 , J) 

DO 21 I=II,IG 
PZ=0 . 0 
PZ1=0.Q 

DO 4018 K= 1 , K MM1 
KP 1 =K +1 

PZ= PZ-ZZ (KP1) *AL2 * (TT (I ,J, K) + TT (I, J,KP1) ) 

PZ1=PZ1 + PZ*Z (KP1 ) 

4018 CONTINUE 

SP(I , J) =P21*HG 
21 CONTINUE 
4011 CONTINUE 
20 CONTINUE 
DO 4019 

SP { I MP1 , J) =SP (3, J) 

4019 CONTINUE 

DO 3000 J=J MH , JMM 1 
F JU J=FJU (J) 

IF (ABS (FJUJ) - LE. 1 -0E-7) GO TO 3 
JP 1 = J -n 
R XM J= R*X M (J) 

R 1 = F XMD* FUU J 
RX=RY*FJUJ 
GO TO 4 

3 JJ=J 

GO TO 3000 

4 IJ=LV(J) , 

DO 4 0 02 M 1 = 1 , I J 

II=LA (Ml, J) 

IG= LB (M 1 , J) 

DO 3001 I«II,IG 
IP1=I+1 

SP0 = SP (1,0) - SP (IP1, JP1) 

SP1 = SP (IP1 , J) “SP (I, DPI) 

PXBIJ=SP 1-SPO 
PYBI Js-SP1-SP0 
DO 3002 K=1,KM 
22K=ZZ(K) 

KM1=K-1 

IF (K - NE. 1) GO TO 4000 
PA-0 . 0 
PB= 0. 0 
GO TO 4001 
4000 CONTINUE 

PI =T (I, JP1 , K) +T ( I , JP1 , K Ml ) -T (IP1,J,K) -T (IP1,J,KM1) 

P2-T{I,J/K)+T {I,J,KM1)~T(IP1,JP1,K)-T (IP1,JP1,KM1) 

Tin n’ •77V. , -»T-. { d 1 n->\ 



4001 


3002 

3001 

4002 

3000 


6 

7 


4003 

4007 


700 1 


7000 

7003 


7 02 4 

^p)2 


7010 


7012 
701 1 


CONTINUE 

PV= ( PA* AQ-PXBIJ) /Hi 
P0= {-PB*AQ+PY BIJ) /RX 
U (I, J,K) = PU 
DU ( 1 , 0, K) = PU 

V {I, 0,K) = PV 
VV(I,U,K)=PY 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 

IF(JJ.EQ.O) GO TO 7 

JM1=JJ-1 

UP1=JJ+ 1 

DO 5 1= I M H , T M 

IF (KQQ (I, JJ) . EQ.2) GO TO 5 
DO 6 K=1,KM 
PA=UU (I, JK1,K) 

PB= UU (I, JP1,K) 

PU=VV (I , JM1,K) 

PV=VV (I,JP1,K) 

IF (PA.GT.200. 0) PA=0*0 
IF(PB.GT.200. 0) PB = 0. 0 
IF (PU-GT.200. 0) PU=G.O 
IF (PV.GT. 200. 0) PV=0.0 
AU= (PA+PB) /2 . 0 
A V- (PU+PV) /2 . 0 
UU (I, JJ,K)=AU 
0(1, JJ,K) =AU .* 

VV (I , JJ # K) = AV 
V(I,JJ,K)=AV 
CONTINUE 
CONTINUE 
DO 4007 K= 1, K M 
DO 40 08 0=0 m , J MM 1 
U {2, J,K) =U (IM,J,K) 

V (2 , 0 , K) =V (IM,J,K) 

UU { 2 , J, K) =UU (IM, J,K) 

VV (2 ,J, K) = VV (IN, J, K> 

U (IMP 1, J, K)=0 (3, J , K) 

V (IMP1, J,K)=V (3, J , K) 

UU ( I MP1 , J, K) =UU (3 , J, K) 

VV{IMP1, J,K)=VV (3 , J, K) 

CONTINUE 

CONTINUE 
CALI CONTI 
DO 7002 K=2,KH 
WHITE (6,7001) K 

FORM AT (1H ,23 HVERTICAL VELOCITY AT K=I1) 
DO 7003 X- INH ,1 M 

WRITE (6,7 000) I, (W (I, J,K) ,J=3,2 0) 

FORM AT ( 1H ,13 , IX, 18F7.2) 

CONTINUE 

DO 7024 I=IMH,IM 

WRITE (6,7000) I, (W(I,J,K) ,J=21, 38) 

CONTINUE 

CONTINUE 

DO 7009 K=1,KM 

WRITE (6,7010) K 

FORM AT ( 1 H ,7HU AT K=I1) 

DO 7011 I=IKH,IM 

WRITE (6,7012) I, (UU(I,J,K) ,J=3,20) 

FORM AT ( 1H ,I3,2X, 18F5.0) 

CONTINUE 

DO 7022 I=IMH,IM 


Si 



WRITE (6,7 012) I , (UU {I , J,K) , J=21 ,3 8) 
7Q22 CONTINUE 

1)0 7021 I=IHH,I» 

* WRITE (6,7 012) 1 , (VV{I,J,K) ,J=3,20) 
7021 CONTINUE 
DO 7023 

WRITE (6,7012) I, {VV (I , J,K) ,J=21,3 8) 
7023 CONTINUE 
7009 CONTINUE 
BETUIN 
END 



